
for the compressible Navier-Stokes and viscous MHD equations 



For a general class of hyperbolic-parabolic systems including the compressible Navier- 
Stokes and compressible MHD equations, we prove existence and stability of noncharac- 
teristic viscous boundary layers for a variety of boundary conditions including classical 
Navier-Stokes boundary conditions. Our first main result, using the abstract framework 
established by the authors in the companion work ^GMWZ6J , is to show that existence 
and stability of arbitrary amplitude exact boundary-layer solutions follow from a uni- 
form spectral stability condition on layer profiles that is expressible in terms of an Evans 
function (uniform Evans stability). Whenever this condition holds we give a rigorous 
description of the small viscosity limit as the solution of a hyperbolic problem with 
"residual" boundary conditions. Our second is to show that uniform Evans stability for 
small-amplitude layers is equivalent to Evans stability of the limiting constant layer, 
which in turn can be checked by a linear-algebraic computation. Finally, for a class 
of symmetric-dissipative systems including the physical examples mentioned above, we 
carry out energy estimates showing that constant (and thus small-amplitude) layers 
always satisfy uniform Evans stability. This yields existence of small-amplitude multi- 
dimensional boundary layers for the compressible Navier-Stokes and MHD equations. 
For both equations these appear to be the first such results in the compressible case. 
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1 Introduction 



In this paper, we study existence and stability of noncharacteristic viscous boundary 
layers of hyperbolic-parabolic systems of the type arising in fluid and magnetohydrody- 
namics (MHD). In a companion paper [GMWZH] . we have shown under mild structural 
assumptions that for such layers, maximal linearized stability estimates, transversality of 
layer profiles, and satisfaction of the uniform Lopatinski condition by the associated residual 
hyperbolic system all follow from a uniform spectral stability condition on layer profiles that 
is expressible in terms of an Evans function (uniform Evans stability, Definition 13. 4p . Here 
we use these abstract results to obtain existence and stability in interesting physical appli- 
cations. Our structural hypotheses are general enough to allow van der Waals equations of 
state. 

In contrast to our previous joint papers, our main concern here is to give an analytic 
verification of the uniform Evans hypothesis for small amplitude layers (after construct- 
ing such layers), and to give explicit calculations for important physical systems (isentropic 
Navier-Stokes, full Navier-Stokes, and viscous MHD). For the isentropic Navier-Stokes equa- 
tions we are also able to construct transversal large amplitude layers, and to check maximal 
dissipativity of the associated residual hyperbolic boundary conditions. 

The main results of this paper are as follows: (i) assuming the uniform Lopatinski condi- 
tion and transversality of layer profiles (Definitions 11.161 and II. lip , we construct arbitrarily 
high-order approximate boundary-layer solutions matching an inner boundary-layer profile 
to an outer hyperbolic solution; (ii) assuming uniform Evans stability, we use results of 
|GMWZ4l IGMWZ6] to show existence and stability of exact boundary-layer solutions close 
to the approximate solutions, and consequently we obtain convergence of viscous solutions 
to solutions of the residual hyperbolic problem in the small viscosity limit; (iii) we show that 
uniform Evans stability of small-amplitude boundary layers is equivalent to uniform Evans 
stability of the associated limiting constant layer; and (iv), we use (iii) to verify the uniform 
Evans condition for small amplitude layers for a class of symmetric- dissipative systems that 
includes the above physical examples as well as the class introduced by Rousset |R3j . 

In connection with (iii) and (iv) above, we prove existence of small-amplitude layer 
profiles for a variety of boundary conditions, including mixed Dirichlet-Neumann conditions. 
These profiles appear in the leading term of the approximate solutions referred to in (i). 

The above results yield existence and stability of multi-dimensional small-amplitude non- 
characteristic boundary-layer solutions of the compressible Navier-Stokes and viscous MHD 
equations. In both cases these appear to be the first such results. For large-amplitude lay- 
ers, the questions of existence and stability are reduced to verification of the uniform Evans 
condition. Efficient numerical methods for such verification are presented, for example, in 
[CHNZj |HLyZ2l |HLyZ2| . 

Although most of the important physical examples can be written in conservative form, 
the general theory is presented for the more general nonconservative case in sections [21 [3l H] 
and the appendices. As we noted in [GMWZ61 IGMWZ7| . the theory is clearer and in many 
ways simpler in the more general setting. 
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For results on stability of boundary layers for the incompressible Navier-Stokes equations 
we refer to [TWllTS]. 



1.1 Equations and assumptions 

Consider as in |GMWZ6j a quasilinear hyperbolic-parabolic system 

d d 

(1.1) Ce{u) := AQ{u)ut + - e ^ dj{Bjk{u)dku) = 0, 

j=l j,k=l 

on [—T, T] X Q, where C is an open set. We assume the block structure 

(1.2) ^o(^) = (l^^ 4')' ^''^''^=(^0 4l)' 
a corresponding splitting 

(1.3) tx = (ni,u2) eM^-^' xR^', 
and decoupled boundary conditions 

T2iu%^dn = 92{t,x), 
T^{u,dTU^,dyU^)\^(,Q^ = 0, 

where Ot and du denote tangential and inward normal derivatives with respect to 50 and 
T^{u,dTU^ -.duV?) = KydyV? + ^^j=iKj{u)Vju'^^ where the Vj{x) are smooth vector fields 
tangent to and Ky is constant. Unless otherwise noted we take bounded with smooth 
(that is, for k large) boundary, but our results apply with no essential change to other 
situations such as the case where dVt coincides with a half-space outside a compact set. 

Setting e = in (jl.ip we obtain Cq, a first-order operator assumed to hyperbolic. The 
parameter e plays the role of a non-dimensional viscosity and for e > 0, the system is 
assumed to be parabolic or at least partially parabolic. Classical examples are the Navier- 
Stokes equations of gas dynamics and the equations of magneto- hydrodynamics (MHD). 

We set 

(1.5) Aj = Aq Aj, Bjii = Aq Bjk, 



d d 

(1.6) A{u,0 = J2^^M^^ ^"^d :B(^x,e) = 5^ CiCfc^ifc(n), 

j=l j,k=l 

and systematically use the notation M"^ for the sub-blocks of a matrix M corresponding 
to the splitting u = (u^,u^). Note that 

(1.7) B,^k{u) := Ao{u)-'B^k{u) = ( -22, . ] , 

\0 Bjkiu), 
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so it is natural to define the high-frequency principal part of (jl.ip by 



(li 



dtu^ + A^^{u,d)u^ = 0, 



Our main structural assumptions are modeled on the fundamental example of the Navier- 
Stokes equations with general, possibly van der Waals type equation of state. For applica- 
tions it is important to allow the states assumed by solutions of (jl.ip when e > or when 
e = to vary in overlapping but not necessarily identical regions of state space. More- 
over, these regions must be allowed to depend on {t,x). These considerations motivate the 
following definition of U, Uq, and U* . 

For some T > let 0{t, x) be a continuous set-valued function from [— T, T] x to open 
sets in M^, and define graphs 

U = {(t, X, 0{t, x)) : (t, x) G [-T, T] x ^] 

Ua = {{t,xo,0{t,xo)) : {t,XQ) e[-T,T]x d^]. 

For (t,xo) G [— r, T] X dVL let 0*{t,XQ) be another continuous op en-set- valued function 
satisfying O*{t,xo) D 0{t,XQ) and define 

(1.10) U* = {{t, xo, 0*it, xo)) : (t, xo) G [-T, T] x dQ}. 
Observe that we have 

(1.11) UaCUnW, 

but neither U nor U* is a subset of the other. For elements of U define 

(1.12) Tr{t,x,0{t,x)) =0{t,x) 

and define tt similarly for elements Ug and U*. Finally, denote by ttU (resp. irU*, ttUq) 
the union of the open sets obtained by applying vr to elements of U (resp. U* , Uq). 

The set t^U is the "hyperbolic set" where solutions of the inviscid equation Cq{u) = 
take their values; t:IA* is the set where boundary layer solutions u'^ of the viscous equations, 
restricted to a small neighborhood of the boundary, take their values. In particular, the 
layer profiles (Definition 11.141) take values in ttU* . The set ttIAq C ttIA H ttU* is the set of 
profile endstates where matching of the two types of solutions occurs; more precisely, it 
contains the limits as 2; ^ oo of layer profiles (see (jl.l5p ). or equivalently, the boundary 
values of solutions of the associated residual hyperbolic problem (see I1.27P . 

Assumptions 1.1. 

(HI) The matrices Aj and Bjk are C°° N x N real matrix-valued funtions of the variable 
u e ttU U ttU* C R'^ . Moreover, for all ueirUUTrU*, det Ao{u) / 0. 

(H2) There is c > such that for all u G ttU U ttU* and G M*^, the eigenvalues of 
B (m,0 satisfy Re/x > c|^p. 
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(H3) For all u £ irU U irW and ^ £ M'^\{0}, the eigenvalues of A (u,^) are real, 
semi-simple and of constant multiplicity. Moreover, for u S -kU* , det ^^'^(n, i/) 7^ 0, with 
the eigenvalues of the normal matrix A (u, i^) 7^ all positive (inflow) or all negative 
(outflow), where v denotes the inward normal to dO,. 

(H4) For allu G ttU and S M'^\{0}, the eigenvalues o/A(n,^) are real, semisimple, and 
of constant multiplicity. Moreover, for u G ttUq, det A{u,i') 7^ 0, with number of positive 
(negative) eigenvalues of A{u,i') independent of u. 

(H5) There is c> such that for u G ttU and G W^, the eigenvalues ofiA{u,^)+B{u,^) 
satisfy Re /x > c ^^'^p . 

Remark 1.2. l.)In Hypothesis (H4) the statement "for u G ttUq, det A^Ujv) ^ 0" should 
be interpreted as asserting that for {t,xo,u) G Uq, we have det A{u,i'{xo)) 0. A similar 
remark applies to (H3) and to later statements of this sort. 

2. ) Hypothesis (H4) is a hyperbolicity condition on the inviscid equation Cq{u) = 0, 
while (H2),(H4) implies hyperbolic-parabolicity of the viscous equation Cs{u) = when 
e > 0. (H3) is a hyperbolicity condition on the first equation in (ll.Sp . The conditions on 
the normal matrices in (H3)-(H4) mean that the boundary is noncharacteristic for both the 
inviscid and the viscous equations. Hypothesis (H5) is a dissipativity condition reflecting 
genuine coupling of hyperbolic and parabolic parts for u G ttU. 

3. )Later we will occasionally drop the vr on irU in statements like u below. 

Symmetry plays an important role in applications such as those to the Navier-Stokes 
and MHD equations considered here. In particular, (H5) holds always when the conditions 
in the following two definitions are satisfied [KaSlj IKaS2j . 

Definition 1.3. The system (ll.ip is said to be symmetric dissipative if there exists a 
real matrix S{u), which depends smoothly on u £ nU, such that for all u G nU and 
all C e K'^\{0}, the matrix S{u)Ao (u) is symmetric definite positive and block-diagonal, 
S{u)A{u,^) is symmetric, and the symmetric matrix Ke S{u)B{u,S^) is nonnegative with 
kernel of dimension N — N' . 

Given that A^ and the matrices Bj^ have the structure (jl.2p . observe that we have the 
equivalences: 

SAq block diagonal <^ S lower block triangular <^ SB block diagonal^ 

Definition 1.4. A symmetric-dissipative system satisfies the genuine coupling condition 
if for all u G ttU and all ^ G M'^\{0}, no eigenvector of ^jCj ^^^^ the kernel of 

The constant multiplicity condition in Hypothesis (H4) holds for the compressible Navier- 
Stokes equations whenever A{u, ^) is hyperbolic. We are able to treat symmetric-dissipative 
systems like the equations of viscous MHD, for which the constant multiplicity condition 
fails, under the following relaxed hypothesis. 

^ The block-diagonal assumption repairs a minor omission in |GMWZ6| : this is needed to conclude (H5) 
from symmetric dissipativity plus the genuine coupling condition. 



6 



Hypothesis H4'. For all u € ttU and ^ G M*^ \ {0}, the eigenvalues of A{u,^) are 
real and are either semisimple and of constant multiplicity or are totally nonglancing in the 
sense of |GMWZ6] . Definition 4.3. Moreover, for u G ttUq we have det A{u, v) ^ 0, with the 
number of positive (negative) eigenvalues of A(u, v) independent of u. 

Remark 1.5. The condition of constant multiplicity in (H3) can probably be dropped for 
symmetric-dissipative systems. For sufficiently small- amplitude boundary layers, we expect 
that the condition in (H3) that eigenvalues have a common sign may be dropped as well; see 
Remark \4-S\ 2. 

Notations 1.6. With assumptions as above, (constant) denotes the number of positive 
eigenvalues of Ai,{u) := A{u,v) for u G ttUq and Nl_ the number of positive eigenvalues of 

aI\u) := A^\u, u) for ueirU*. We also set Nb = N' + Nl. 

As indicated by block structure (|1.8|) . is the correct number of boundary conditions 
for the well posedness of (jl.ip . for solutions with values in U* U U, with A'^' boundary 
conditions for and boundary conditions for u^. On the other hand, iV+ is the correct 
number of boundary conditions for the inviscid equation for solutions with values in U. 

Assumption 1.7. (H6)Ti, T2 and T3 are smooth functions of their arguments with values 
inR^+, R^'~^" and R'^" respectively, where N" G {0, 1, . . . , iV'}. Moreover, has 
maximal rank N" and for allu G U* the Jacobian matrices T'^iu}-) and ^^{u'^) have maximal 
rank N]_ and N' — N" respectively. 

Examples 1.8. Hypotheses (H1)-(H5) are satisfied by the compressible Navier-Stokes 
equations, provided the normal velocity of the fluid is nonvanishing on Li* and the normal 
characteristic speeds (eigenvalues of A{u,v)) are nonvanishing on Uq. We have N]_ = 1 
or according as normal velocity is positive with respect to inward normal (inflow) or 
negative (outflow). Recalling that li* is the set where boundary layer solutions of the 
viscous equations, restricted to a small neighborhood of the boundary, take their values, we 
see that these velocity restrictions correspond to having a porous boundary through which 
fluid is pumped in or out, in contrast to the characteristic, no-flux boundary conditions 
encountered at a solid material interface for which normal velocity is set to zero. 

Hypotheses (H1)-(H5), with (H4) replaced by (H4'), are satisfied by the viscous MHD 
equations with ideal gas equation of state under similar velocity restrictions on the plasma. 
See Remark 15.61 and the discussion following it for the precise velocity requirements and for 
examples of boundary conditions for viscous MHD and the corresponding reduced hyper- 
bolic problem. 

We note that our structural hypotheses are sometimes satisfied under much weaker 
assumptions on the equation of state, which may be of van der Waals type on U* and just 
thermodynamically stable on U ; see |GMWZ4l IZ3| . 

Boundary conditions of the type we consider for the compressible Navier-Stokes equa- 
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tions are of considerable practical importance. Recall the equations are 



' dtp + div(/9u) = 



(1.13) 



< 



dt{pu) + div(/9n*n) + Vp = epAu + e{p + r/)Vdivn 
dtipE) + dw{{pE + p)u) = kAT + e^div((n • V)u) 



+ e{p + T])'V{u ■ divn) 



where p denotes density, u velocity, e specific internal energy, E = e + specific total 
energy, p = p{p, e) pressure, and T = T{p, e) temperature. Take the unknowns to be 
{p, u, T), and consider the problem on the exterior = a'^ of a bounded set a with smooth 
boundary, with no-slip suction-type boundary conditions on the velocity, 

urlan = 0, u^lan = V{x) < 0, 

and either prescribed or insulative boundary conditions on the temperature, 

T\dn = Ty,au{x) or di,T\Q^ = 0. 

This corresponds to the situation of an airfoil with microscopic holes through which gas is 
pumped from the surrounding flow, the microscopic suction imposing a fixed normal veloc- 
ity while the macroscopic surface imposes standard temperature conditions as in flow past 
a (nonporous) plate. This configuration was suggested by Prandtl and tested experimen- 
tally by G.I. Taylor as a means to reduce drag by stabilizing laminar flow; see [HI IBrj . It 
was implemented in the NASA F-16XL experimental aircraft program in the 1990's with 
reported 25% reduction in drag at supersonic speeds [BrjH 

Remark 1.9. At the expense of further bookkeeping, we could equally well define separate 
Uqj C U* for each connected component {di})j of the boundary di}, each with distinct values 
of N^, N^. where now the sets Uq j are to be thought of as the sets of possible boundary 
values for the hyperbolic problem at each {dQ)j . This would allow, for example, the situation 
that fluid is pumped in through one boundary and out another as considered in JTW^ for the 
incompressible case. 

We follow hypotheses (II1)-(II6), in some cases with (H4') in place of (H4), throughout 
the paper; unless otherwise indicated, they are assumed in all statements and propositions 
that follow. 

1.2 Layer profiles, transversality, and C manifolds 

To match solutions of the inviscid problem approaching a constant value u at xq £ dO, 
to solutions satisfying the hyperbolic-parabolic boundary conditions, one looks for exact 
solutions of (jl.ip (II. 4[) on the half-space (x — xq) ■ z^(xo) > tangent to at xq of the form 



(1.14) 




2 



See also NASA site http://www.dfrc.nasa.gov/Gallery/photo/F-16XL2/index.html 
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lyixo) the inward unit normal to dO, at xq, such that 
(1.15) Um w{z) = u. 

2:— >+oo 

Solutions are called layer profiles. 
The equation for w reads 



(1.16) 



A^{w)d,w - d,{B4w)d,w) =0, z > 0, 
T(u;,0,5^?i;^)|^=o = {9iit,xo),g2{t,xo),0), 



where Ai,{u) = X]j=i By{u) := Bjk{u)vjVk- The natural limiting boundary 

conditions for the inviscid problem should be satisfied precisely by those states u that are 
the endstates of layer profiles. This leads us to define 

(1.17) C(t,xo) = {u : there is a layer profile w G C°°{W+;U*) satisfying (fTTS]) . (fTTe]) . 

The profile equation (11.16P can be written as a first order system for U = {w,dzw'^), 
which is nonsingular if and only if A}} is invertible, (H3): 

d.w' = -{Al'r'Al'w^ 

(1.18) d,w'^ = 

d.{Bl'w') = {Af-Al\Al'r'Al')w\ 

and the matrices are evaluated at it; = {'w^,w'^). 

Consider now the linearized equations of (I1.16|) about w^z), written as a first-order 
system 

(1.19) d,W - g^{z)W = 0, z>0, 

(1.20) r,VF|,=o = 
inW = {w^ , w'^ , w'^) , where 

/o -iK'r'Ai^ 

(1.21) Gui+oo) := lim G^iz) =00 / \ (u) 

\0 (i322)-l(^22_^21(^ll)-1^12, 

and (note decoupling between and (n^,n^) := {v^^jdyV?) variables) 

(1.22) TyU =[Tru^,T2U^,K^u^). 



9 



Lemma 1.10 ( |GMWZ"6| IMaZ3] ) . Let N'^ denote the number of stable eigenvalues 3ft/i < 
of Gui+oo), the number of unstable eigenvalues "Sifi > 0, 5 the subspace of solutions of 
(|1.19p that approach finite limits as z ^ oo, and Sq the subspace of solutions of (|1.19p that 
decay to 0. Then, 

(i) Nl + Nl = N' and 

(1.23) N+ + nI = Nb := N' + nI, 

(a) profile w{-) decays exponentially to its limit u as z ^ +00 in all derivatives, and 
(Hi) dim5 = N + Nl and dim^o = Nl. 

Proof. A proof is given in [GMWZH] . Lemma 2.12. □ 

Definition 1.11. The profile w is said to be transversal if 

i) there is no nontrivial solution w £ Sq which satisfies the boundary conditions 
Ti,{w,dzW^)l^=o = 0, 

a) the mapping w 1— > riy{w,dzw'^)\z=o from S to C^*" has rank N),. 

Definition ILllT i) and (ii) corresponds to the geometric conditions that the level set 
{W : T{W) = T(u;(0), 0, 9zW^(0))} have transversal intersections in phase space W = 
{w,dzUp') at Wq := {w{&),dzUp'{0)) with the stable and center-stable manifolds, respec- 
tively, of VFoo := (M)0); see Lemma [53] below. 

The following assumption is the starting point for our construction of exact boundary 
layer solutions to (jl.lh . 

Assumption 1.12. Fix a choice 0/(51,(72) o-s in (II. 4p . ForUg as in ()1.9p we are given a 
smooth manifold C defined as the graph 

(1.24) C = {{t,xo,C{t,xo)) : {t,xo) e [-T,T] x OQ} C U9, 

where eachC{t,X()), defined as in (jl.l7p . is now assumed to be a smooth manifold of dimen- 
sion N — . In addition we are given a smooth function 

(1.25) : [0, 00) X C ^ vrZ^* 

such that for all {t,xo,q) £ C, W{z,t,xo,q) is a transversal layer profile satisfying (I1.16P 
with V = i'(t,xo) and converging to q as z ^ 00 at an exponential rate that can be taken 
uniform on compact subsets ofC. 

This assumption is hard to check in general. However, in Proposition 12.61 we show that 
for a large class of problems including the Navier-Stokes and MHD equations with vari- 
ous boundary conditions, Assumption 1 1 . 1 2) is always satisfied for small-amplitude profiles. 
In Proposition 12.41 we give necessary and sufficient conditions on boundary operators of 
the form (jl.4p in order for Assumption 11.121 to hold in the small-amplitude case. These 
boundary conditions include the standard noncharacteristic boundary conditions for the 
Navier-Stokes and viscous MHD equations. In Proposition 12.81 we give a local construction 
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of a C manifold with associated profiles near a given, possibly large amplitude, transversal 
profile. In Proposition 12.101 we give a global construction of a C manifold with associated 
profiles near a given family of large-amplitude profiles, assuming such a family exists (not 
clear in general). 

As regards the failure of Assumption I1.12[ Proposition 12.61 implies, for example, the 
following nontransversality result. 

Proposition 1.13. Whenever rankTs = N" > N"^, small- amplitude profiles are not 
transversal. Equivalently, whenever the number of scalar Dirichlet conditions for the parabolic 
problem, (N' — N") + N^, is strictly less than the number of scalar boundary conditions for 
the residual hyperbolic problem, (see ()1.27p ). small- amplitude profiles are not transver- 
sal. 

1.3 Inviscid solutions, residual boundary conditions, and the uniform 
Lopatinski condition. 

Under suitable assumptions we will study the small viscosity limit of solutions to 

d d 
Ce{u) := AQ{u)ut + ~ dj{Bjk{u)dku) = 0, 

T{u,dTu\d^u^) = (51,52,0) on [0,To] x dfl 

and demonstrate convergence to a solution u^{t,x) of the inviscid hyperbolic problem: 

Co(u°) = on [0,To] x n 
(1.27) ' ^ ^ 

n°(t, xo) G C{t, xo) for (t, xq) G [0, Tq] x dn, 

where C{t,XQ) is the endstate manifold defined in Assumption 11.12] (see also Prop. 12. 6p . 

Definition 1.14. We refer to the boundary condition in (jl.27p as the residual hyperbolic 
boundary condition. 

For (t, Xq) G [0, Tq] X dQ we freeze a state p := u^{t,xo) and, working in coordinates 
where the boundary is = 0, we define 

(d-l 
{iT + ^)Ao{p) + J2iVjAj{p) 

Here we have suppressed the dependence of the frozen Aj on spatial coordinates in the 
notation. Let 

(1.29) ip:R^ ^ 
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be a defining function for C{t,xo) near p, i.e., C{t,xo) = {u : ip{u) = 0}, with Vip full 
rank A^_|_. Then, the residual boundary condition (jl.27p may be expressed, locally to p, as 
'^resiu) := ip{u), hence the linearized residual boundary condition at p takes the form 

(1.30) Tres{p)u = <^ i;'{p)u = <^ it G TpC{t, xq). 

Remark 1.15. Suppose w{z) is a solution of ()1.16p converging to p = u^{t,XQ) G C(t,xo) 
as z ^ CO. Let us write the linearized equations of (I1.16P around w{z) as 

(1.31) l.{t,xo,z,dz)w = 0, T^(^^^){w,w'i) = 0. 

Observe that the tangent space TpC(t, xq) niay be characterized as the set of limits at z = oo 
of solutions to (|1.3ip . This follows readily from the definition of C{t,xo); see jMet^, Prop. 
5.5.5. 

Definition 1.16. 1) The inviscid problem (|1.27p satisfies the uniform Lopatinski condition 
at p = u{t, Xq) provided there exists C > such that for all (" with 7 > 

(1.32) \DLop{p,C)\ := |det (E-(F(p, C)), ker r,e.(p)) | > C. 

2) The inviscid problem (jl.27p satisfies the uniform Lopatinski condition provided (|1.32p 
holds with a constant that can be chosen independently of{t,xo) G [0,To] x dQ. 

Here by a determinant of subspaces we mean the determinant of the matrix with sub- 
spaces replaced by smoothly chosen bases of column vectors, specifying Dlop up to a smooth 
nonvanishing factor; for example, if the bases are taken to be orthonormal, then |-DlopI is 
independent of the choice of bases. We refer to |GMWZ7] . section 4.1 for equivalent formu- 
lations and further discussion of the uniform Lopatinski condition. 

Theorem 1.17. Given a smooth manifold C as in Assumption M.lB, consider the hyperbolic 
problem ()1.27p 

(i) under hypotheses (H1)-(H5), or alternatively, 

(a) assuming (H1)-(H5), except that (H4) is replaced by (H4') in the symmetric-dissipative 
case. 

Let s > I + 1 and suppose that we are given initial data v^{x) G H^^^{i}) at t = 
satisfying corner compatibility conditions to order s — 1 for (|1.27p . Suppose also that the 
uniform Lopatinski condition is satisfied at all points xq G di}, t = 0. Then there exists a 
To > and a function u^{t,x) G H^{[0,Tq] x Q) satisfying (jl.27p with 

(1.33) nf ^0 = v^, 

and so that the uniform Lopatinski condition holds on [0,To] x dQ. 

Proof. We refer to \CP\ . Chapter 7 for a discussion of corner compatibility conditions. For 
the proof in the case of constant multiplicity (i.e., when (H4) holds) see |CPj . Chapter 7 
and |Met2] . For the case of variable multiplicities (including situations more general than 
those considered here) see |MZ2| . 

□ 
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1.4 Approximate solutions to the viscous problem 



The first step in constructing exact solutions to the viscous problem (|1.26|) that converge 
to a given solution of the hyperbolic problem (jl.27p in the small viscosity limit is to 
construct high order approximate solutions of (jl.26p with that property. Following the ap- 
proach of |GMWZ4l [GMWZ7j we construct approximate solutions using a WKB expansion 
as described in the following result. A more precise statement of Proposition 11.181 and the 
proof are given in Appendix Rl 

Let 



(1.34) {xo,z) -.Q^ dnxR+ 

be a smooth map defined for d{x,dn) <r,r>0 sufficiently small, such that 

{xq,z)~'^{xo,z) = Xo + zi^{xo), 

where 1^(2:0) is the inward normal to dO, at xq, i.e., {xq,z) are normal coordinates and 
Vz = on Let xi^) be a smooth cutoff function vanishing for d{x,dQ) > 2r and 
identically one for d{x,dQ) < r. 

Proposition 1.18 (Approximate solutions). Suppose we are given a C manifold and asso- 
ciated transversal profiles W{z, t, xq, q) as in Assumption \l.l^ and also a solution S H'^^ 
to the inviscid problem (I1.27P as described in Theorem (In particular, the uniform 

Lopatinski condition is satisfied on [0,Tq] x dil.). Fix positive integers M and s with 
M > s > 1. Provided sq is sufficiently large relative to M and s (see (IA.25|) ). there 
exists an approximate solution u^{t,x) to the viscous problem (|1.26p on [0,Tq] x Q of the 
form 

ul{t,x)= em\t,x,-) + e^^+^u^"{t,x), 

(1.35) o<J<^^ 

W{t,x,^) = x{x) {wH-^,t,XQ) -W\+oo,t,xo)) +u^{t,x). 

Here is the given inviscid solution, the leading profile is given by 

(1.36) W'^iZ, t, xo) := W{Z, t, xq, u\t, xo)) 

for W{Z,t,xo,q) as in Assumption \l.l^ and the higher profiles {Z,t, xq) converge ex- 
ponentially to their limits as Z ^ -\-oo. The approximate solution satisfies 

\\^e{Ua)\\H''i[0,TQ]xn) < Cs^^ 

T{ua,dTul,duul) = (51,52,0) on [0,To] x 
for (91, (72)0) as in the original viscous problem (|1.26p . 
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Examples 1.19. The combination of Theorem II. 171 (for inviscid solutions u^), Proposition 
12.61 (for C manifolds and transversal profiles), and Corollary 11.291 (for satisfaction of the 
uniform Evans condition which implies transversality and uniform Lopatinski) provides us 
with a variety of examples, involving both Dirichlet and Neumann conditions for the Navier 
Stokes and viscous MHD equations and including all cases mentioned in Examples II. 8^ 
where the hypotheses of Proposition 11.181 are satisfied for small-amplitude profiles. 

The large-amplitude case is more difficult. In section [5711 we construct C manifolds with 
associated large- amplitude profiles for the isentropic Navier-Stokes equations. Transversal- 
ity of those profiles is verified in Proposition 15.51 and maximal dissipativity of the residual 
hyperbolic problem, which implies the uniform Lopatinski condition, is verified in section 
15.1.11 Thus, here again the hypotheses of Proposition 11.181 are satisfied. 

Large-amplitude C manifolds and transversal profiles are constructed locally near the 
endstate of a given transversal profile in Proposition 12.81 In problems with rotational sym- 
metry it may be possible to promote locally constructed large-amplitude C manifolds to 
global ones by using that symmetry. For large-amplitude layers there are now efficient nu- 
merical methods available for verification of the uniform Evans condition. [CHNZl HLyZ2 
HLyZ2| 



1.5 The Evans condition, stability and transversality, and the small vis- 
cosity limit 

We refer to the approximate solution described in Proposition II . 181 as an approximate 
boundary-layer solution. Under a suitable Evans condition (Definition 1 1 . 20p . we will produce 
a nearby exact boundary-layer solution of <\l.lh . (II. 4p of the form 

(1.38) = ul + v', 
where satisfies the "error problem" 

(1.39) Ce{Ua +v) - Ce{Ua) = -CeiUa), 

derived from the problems satisfied by and n^, with boundary data 

(1.40) {Tiua + V, OtA^I + v^)) - T{ua, dT,uul)) Uan = 

and forcing term —Ci;{ua) small of order e^^ in as described in (ll.STh . with M and s 
large. Here t'|t=o satisfies high order corner compatibility conditions depending on Ua- 

Evidently, the problem of estimating convergence error ||v^||_f/s in terms of truncation 
error ([L39]) amounts to determining ff" H'^ stability estimates for (11.39p - (ll.40p . for 
which the main obstacle is rapid variation of coefficients in the boundary-layer region. We 
focus now on this region, and stability of associated layer-profiles, i.e., the "microscopic" 
stability problem. 

Fix a point {t, xq), xq G d^l and consider again the viscous problem (11.26p . We work in 
local spatial coordinates {y, Xd) where xq is (0, 0) and dO, is given by Xd = 0. Consider a 
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planar layer profile 

(1.41) u''{t,y,Xd) =w{xd/e) 

as in (jl.l4p . which is an exact solution to (jl.26p on > when the coefficients and 
boundary data (51, 52)0) are frozen at (t, xq). Without loss of generality we take £ = 1, set 
z = Xd, and write the linearized equations of (jl.26p about w: 

(1.42) C'^u = f, T'{u,dyu\d,u%=o = 9- 

Here T' is the differential of T at {w{0),0,dzw'^{0)) and £^ is a differential operator with 
coefficients that are smooth functions of x^- 

Performing a Laplace- Fourier transform of (jl.42p in (t, y), with frequency variables 
denoted by 7 + ir and ry respectively, yields the family of ordinary differential systems 

(1.43) L{z,-/ + iT,ir],dz)u = f, T' {u,ir]u^ ,dzU^)\z=Q = g, 

(1.44) L = -B{z)dl + A{z, Qd, + M{z, C), 

with in particular, B{z) = Bd4{w{z)) and A^^{z,C,) = -^]i{w{z)). The matrices A{z,C,): 
M.{zX) are written out explicitly in ()B.30p . 

The problem (|1.43p may be written as a first order system 

(1.45) d,U -G{z,Q)U = F, r(C)[/|,=o = G, 

where U = {u,dzu'^) = (u^jU^jU^) G C^+^' and C = iT,^,rj). The components of G{zX) 
are given explicitly in (IB.33P and we have 

(1.46) T{C)U = (r^ni, r^n^, r3(C)n3) 
with 

r'u' = T[iw\0))n\ rV = r,{w\o))u\ 
T^{ri)U = KdU^ + KT{w{0))iriu^. 

Observe that when C = (I1.45P coincides with (11.19| )- (jl.20[) in the case when i' = (0, 1). 

A necessary condition for stability of the inhomogeneous equations (jl.45p is stability of 
the homogeneous case F = 0, G = 0, i.e., nonexistence for 7 > 0, 7^ of solutions U 
decaying as 2; — > +00 and satisfying r(C)C/(0) = 0. These may be detected by vanishing of 

the Evans function 

(1.48) D{C) := ^d_et/E-(C),kerr(C)), 

where lE^(C) is the subspace of initial data at z = for which the solution of dzU—Q{z, C,)U = 
decays at z = +00. For high frequencies |C| > i? > we also define in (|3.14p a rescaled 
Evans function Z)*'^(("). 

Given a C manifold and associated layer profiles W{Z,t,XQ,q) as in Assumption 11.121 
along with an inviscid solution as in Theorem II. 17| we define in the same way D{t, xq, C) 
and D'"^(t, xq, C) for every {t,xo) £ [0, Tq] x 917 using the associated profile ti;(Z) = W{Z,t,XQ,u{t,XQ)). 
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Definition 1.20. 1.) We say that the uniform Evans condition is satisfied by the profile 

w{z) provided there exist positive constants R, C such that 



1^(01 >C /or 0< Id <i?, 7>0 and 
^ ■ ' \D'\0\>C for\C\>R,i>Q. 

2.) We say that the uniform Evans condition holds on [0, Tq] x provided there exist pos- 
itive constants R, C such that the estimates (jl.49p hold for the Evans functions D{t,xo,() 
and D^^{t,XQ,C) defined above, uniformly for all (t, xq) G [0, Tq] x dQ. 

In order to understand the behavior of D(C) near ^ = it is helpful to introduce polar 
coordinates 

(1.50) ( = pC, /9=|C|, forC/0,7>0 
and to write 

(1.51) D{C) = D{C,p), E-(C)=E-(C,p) forC/0,7>0. 

It is shown in jMZ3j . Theorem 3.3 and [GMWZ6] . Remark 2.31 that the spaces E^((^) = 
E~((^, p) have continuous extensions to /) = under our structural hypotheses (and, more 
generally, when there exist i^— families of symmetrizers for the linearized viscous problem.) 
Hence, under our assumptions on F, D{(,p) extends continuously to p = for ^ with 
7 > 0. This continuity allows us to rephrase the low frequency uniform Evans condition 
equivalently as 

(1.52) D(C,0)/Ofor |C| = 1, 7>0. 

The following elementary result allows us to verify Evans conditions by proving trace 
estimates. 

Lemma 1.21 ( |GMWZ6| . Lemma 2.19). Suppose that E C C" and L : C" ^ C™, with 
rank T = dimE = m. If \ det(E, kerF)! > c > 0, then there is C, which depends only on c 
and |F*(FF*)-i| such that 

(1.53) \U\ < C\rU\ for all U £ E. 

Conversely, if this estimate is satisfied then | det(E, ker F)| > c > 0, where c > depends 
only on C and |F|. 

Remark 1.22. By Lemma \1.21\ the uniform Evans condition |-D(C)| > C > on some 
subset S of frequencies is equivalent to 

(1.54) 1^71 < C|F;7| for ah U G E^(C) 
for some constant C > independent of ( G S. 
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We now recall two results from [GMWZH] . The first extends an earlier result of Rousset 
[R2) and shows that only the low frequency Evans condition is needed for the construction 
of high order approximate viscous solutions w^. 

Lemma 1.23 ([ GMWZ6j Theorem. 2.28). Assume (H1)-(H6) (with (H4') replacing (H4) 
in the symmetric- dissipative case), and consider a layer profile w{z) — > p as z ^ oo. The 
uniform Evans condition holds for low frequencies, that is, there exist positive constants r, 
c such that 

(1.55) \D{C)\ > c for \C\ < r, 

if and only if w is transversal and the uniform Lopatinski condition (Definition M.lb]) holds 
at p for the residual hyperbolic problem (jl.28p - (jl.30p . 

The next result shows that the full uniform Evans condition implies maximal linearized 
stability estimates for the viscous problem. 

Proposition 1.24 ( |GMWZ6j . Theorems 3.9 and 7.2). Assume (H1)-(H6) (with (H4') 
replacing (H4) in the symmetric dissipative case) and consider the problem (jl.42p defined 
by linearization around a layer profile w{z). If the uniform Evans condition is satisfied 
by w{z) f Definition \1.20\) . then solutions to (jl.42p satisfy the maximal stability estimates 
(ffrH]) and (13:81) . 

Proposition 11.241 is proved by constructing smooth Kreiss-type symmetrizers for the 
Laplace-Fourier transformed equations (jl.45p in the low-, medium-, and high-frequency 
regimes. This method of proof is just as important as the result, since the same sym- 
metrizer construction may be used to obtain maximal stability estimates for the linearized 
equations about general (nonplanar) solutions of the form ()1.35p . The procedure, which 
is used in [MZH IGMWZ31 IGMWZ4] , is to freeze slow variables in the original linearized 
viscous problem, take the Laplace-Fourier transform to obtain ODEs depending on fre- 
quency like (jl.45p . construct symmetrizers for the transformed problem, and then quantize 
those symbols to produce paradifferential operator symmetrizers for the original (unfrozen) 
problem. The resulting linear estimates can be used to prove convergence of a nonlinear 
iteration scheme that yields existence of an exact solution to the viscous problem ()1.26p 
that is close to the approximate solution u^, as stated in the following theorem. 

Collecting these observations and combining with Lemma 11.231 and Theorem 11.181 we 
obtain the following main result, which reduces the problem of proving existence and non- 
linear stability of boundary-layer solutions to verification of the uniform Evans condition. 

Theorem 1.25. Consider the viscous problem (ll.26|) under assumptions (H1)-(H6) (or 
with (H4' ) replacing (H4) in the symmetric-dissipative case). Given an inviscid solution 
g i?*o([0,To] X Q) as in Theorem suppose that the uniform Evans condition holds 

on [0, To] X do. (Defn. \1.20\) . Suppose the constants k, M, and sq satisfy 

(1.56) A;>^ + 4, M>k + 2, sq > + ^ + 2M + 
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Then there exists Eq > 0, an approximate solution as in (jl.35p satisfying (|1.37p . and an 
exact solution of (jl.26p such that for < e < eo 

II ^ ^^M—k 

If ~ ^alwi'°°{[0,To]xn) ^ ^£ ) 

(1-57) lk-^^°llL2(f,x[o,T])<Cei/', 

- n° = 0(e) in L^^([0,ro] x !^°) 

where Q° denotes the interior of Q. Moreover, the linearized equations about either u^^ or 
satisfy maximal stability estimates. 

Proof. The proof is by the same iteration scheme as the corresponding result for shocks, 
Theorem 6.18 in |GMWZ4] . except that it is simpler because the function tjj defining the 
free transmission boundary in |GMWZ4] is absent in the present fixed boundary problem. 
In the partially parabolic case the iteration scheme, which is explained in |GMWZ4'j section 
6.1.1, must be designed carefully, because the linearized estimates give weaker control over 
the "hyperbolic component" than the "parabolic component" u^. Higher W^'°° norms 
of — can be estimated by increasing k, M, and sq. 

The proof of Theorem 6.18 in [GMWZ4j used the constant multiplicity assumption (H4) 
to construct symmetrizers. The proof of our Proposition 11.241 yields symmetrizers under 
the weaker assumption (H4') when (H4) fails in the symmetric-dissipative case. Those 
symmetrizers are then used exactly as in the constant multiplicity case. 

□ 



1.6 Verification of the Evans condition for small-amplitude layers 

For large-amplitude boundary-layers, the Evans condition may be checked numeri- 
cally; see [CHNZl HLyZl| in the one-dimensional case, |HLyZ2| in the multi-dimensional 



shock case. For small-amplitude layers, it may be verified analytically in several inter- 
esting cases. In particular, we will show that the uniform Evans condition always holds 
for small-amplitude layers of symmetric-dissipative systems with Dirichlet boundary con- 
ditions, and never (by Proposition I2.5f iil together with Lemma ll.23p for constant layers of 
systems with mixed Dirichlet-Neumann conditions for which the number of incoming 
hyperbolic characteristics on U exceeds the number of Dirichlet conditions imposed on U* , 
or, equivalently, the number of Neumann conditions exceeds Ni, — N-\- = N^. 

Our main spectral stability result is the following Theorem established in Section [3l 
which implies that Evans stability of small-amplitude layers w{z) is equivalent to Evans 
stability of the constant-layer limit w = w{oo), a linear- algebraic condition that can in 
principle be computed by hand. This is in sharp contrast to the shock wave case, for which 
the small-amplitude limit is a complicated singular-perturbation problem [Met 31 \PZ\ \FS1\ 
IFS2] . 

Definition 1.26 (Small amplitude profiles). WithUg as in ()1.9p define 
(1.58) Uq^i, = {{u, u{t, xo)) : {t, Xo,u) G Uq}. 
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For e > and any compact set D C Uq^u, the set of e-amplitude profiles associated to D is 
the set of functions w{z) = w{z,u,v) for which there exist (n, z^) E D such that: 

a) Ay{w)dzW — dz{Bi,{w)dzw) = Q on z>Q, 

h) w{z, u, u) ^ u as z ^ oo, 

c) \\iw,w'l) - (u, 0)11^00(0,00) <£,W-k\< £• 
When e is small we refer to such profiles as small amplitude profiles. 

Remark 1.27. Observe that small amplitude profiles are defined without specifying any 
boundary condition at z = 0. We define the Evans function for such a w{z, u, v) using the 
same formula as before (11.481) . where ^~{C} T{C,) are now defined using linearization 
of (jl.26p around w{z,u,iy). 

Theorem 1.28. For any compact subset D C Ug,^ there exists an e > such that the 
uniform Evans condition is satisfied for the set of e-amplitude profiles associated to D 
( Definition 1 1 . 26\) if and only if it is satisfied for the set of constant layers {w(z, u,v) : w = 
u for all z, (u, v) € D}. 

For any bounded set of frequencies the equivalence in the above theorem can be checked 
by showing that the Evans function depends continuously on parameters, including param- 
eters that measure deviation of profiles from constant layers. Thus, the main difficulty is to 
deal with the unbounded high frequency region. Here our strategy, carried out in Theorem 
13.61 is to show that high frequency Evans stability for ()1.45p (the second condition in ()1.49p ) 
is equivalent to high frequency Evans stability of the quasihomogeneous, decoupled, frozen 
coefficient system (j3.17p . By virtue of the quasihomogeneity and frozen coefficients, the 
latter condition can be checked by relatively simple energy estimates for a compact set of 
frequencies, namely the parabolic unit sphere. 

As a corollary we obtain the following result, established by energy estimates in Sec- 
tion m which establishes uniform Evans stability for small amplitude layers in a variety of 
situations. 

Corollary 1.29. (a) In the strictly parabolic case (N = N' ) the uniform Evans condi- 
tion is satisfied for sufficiently small- amplitude layers (in the sense of Definition \1.2B\) of 
symmetric-dissipative systems with pure Dirichlet boundary conditions, 

rankT^ = 0, 

or with pure Neumann boundary conditions when rankT'^ = N = N'^ . 

(b) In the partially parabolic case (N' < N), the uniform Evans condition is satisfied 
for sufficiently small- amplitude layers of symmetric-dissipative systems with pure Dirichlet 
boundary conditions or with mixed boundary conditions when rankT^ = N' = and 

is totally outgoing (N]^ = Oj. 

(c) In the partially parabolic case when rankT^ = N' = and A^y is totally incoming 
(N}^ = N — N'), Evans stability fails in general even for small amplitude profiles (see 
Example \4-3^ . 

(d) The uniform Evans condition fails for sufficiently small amplitude solutions with 
mixed boundary conditions whenever rankT^ > N"^ ; see Corollary \2.5[ 
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Combining Theorem 11.251 Proposition 12.61 and Corollary ll.29| we obtain the follow- 
ing result asserting existence and stability of small-amplitude boundary-layer solutions for 
symmetric-dissipative systems with various types of boundary conditions. Suppose we are 
given a smooth, global assignment of states 

(1.59) {t,xo,p{t,xo)) G Ud for all {t,xo) G [-T,T] x dQ 
satisfying the viscous boundary condition (jl.4jl : 

(1.60) (Ti(pi(t,xo)),T2(p2(t,xo),T3(p(t,xo),0,0)) = {gi{t,xo),g2{t,xo),0) 
for each {t,XQ). Note that the image of p is compact by compactness of dQ. 

Theorem 1.30. Consider a symmetric- dissipative system (11.260 with boundary conditions 
of the type described in Corollary \1.29\ (a),(b) under hypotheses (H1)-(H6), but with (H4' ) 
in place of (H4). Given a smooth global assignment of states p{t,xo) as in (11.590 -( fl.GOp . 
there exists aC manifold satisfying Assumption (jl.l2p with p{t, xq) S C(t,xo) C ttUq for all 
(tjXo), and associated small amplitude profiles W{z,t,XQ,q) satisfying the uniform Evans 
condition on [—T,T] x dQ. The manifold C defines a residual hyperbolic boundary condition 
as in (fOTD . 

Civen initial data satisfying appropriate corner compatibility conditions for the hyper- 
bolic problem (11.271) . there exists an inviscid solution as in Theorem \l.ll\ an approximate 
solution ti^ as in Proposition I i.gJl and an exact boundary layer solution satisfying all 
the conclusions of Theorem \1.25\ for constants sq, k, M as described there. 

1.7 Application to identifying small viscosity limits. 

The exact viscous solutions in Theorems 11.251 and 11.301 are chosen to satisfy high-order 
corner compatibility conditions at t = that depend on the approximate solution it^; the 
construction of depends in turn on having a C manifold with associated layer profiles 
and an inviscid solution to start with. In this section we show how Theorem 11.251 can 
sometimes be used together with Corollary 11.291 and our results on C manifolds in section [2] 
to identify small viscosity limits of solutions to viscous boundary problems like (|1.26p . even 
when neither the C manifold nor the inviscid solution is given in advance. 

Consider the viscous problem ()1.26p on a half-space $7 = {x E M'^ : > 0}: 

d d 
Ce{u) := AQ{u)ut + '^Aj{u)djU - e ^ dj{Bjk{u)dku) = 0, 

(1.61) ^'''=^ 
Tiu,dTu\d,u^) = (51,52,0) on [-T,T] x dn, 

u = u£ ttUq in t < 0, 

where u is a constant state and the (nonconstant) boundary data {g^{t,xo), g'^{t,xo),0) is 
C°°, equal in t < to the constant 

(1.62) T(n, dru, d,u) = T(n, 0, 0) := (g, 0), 
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and also equal to {g, 0) outside a compact set in [— T, T] x Q. We suppose that (jl.26p is 
a symmetric-dissipative system satisfying (H1)-(H6) with (H4') in place of (H4), and that 
the boundary conditions are of the type described in Corollary II. 29r a)(b). 

Set = (0, 1) and note that the constant profile w{z) = u satisfies the uniform Evans 
condition by Corollary ll.29t so in particular, it is transversal. Now apply Proposition 12.81 
to find neiff hborhoods O C of and O C M^*-"^" of g, smoothly varying manifolds 
Cu,g C O for g ^ O, and transversal small-amplitude profiles 

(1.63) W^^,g{-, q) : [0, oo) x C^^g ^ ttU, 
such that for each q G Ci^^g, w = W^^g{-, q) satisfies 

Ay{w)dzW — dz{Bj^w)dzw) = on z > 

(1.64) T(w;,0,w;,2)(0) = (g,0) 

w{z) q as z ^ oo. 

Provided Tq is sufficiently small, our assumptions on the boundary data imply that 

(1.65) {g\t,xo),g\t,xo)) G O for all (t,xo) G [-To, To] x OQ. 
If we now define 

^^gg^ C(f,a;o) :=C^,g(t,^o) and 

W{z,t,xo,q) := W^^g(^t,xo)iz,'l) for {t,xo) G [-Tq^Tq] x dQ, 

the manifolds C(t, xq) and profiles W{z,t,xo,q) satisfy the conditions of Assumption 11.121 
Moreover, by Corollarv 11.291 the uniform Evans condition holds on [— To,ro] x dO,. 

Next, recalling Lemma 11.231 we apply Theorem 11.171 to construct the unique solution 
u^{t,x) G H'^°{[—To,Tq] X rj) to the inviscid hyperbolic problem 

£o(n°) = on [-To,ro] x Q 

(1.67) M°(t, xo) G C{t, xq) on [-Tq, Tq] x dn, 

M° = u in t < 0. 

on a possibly shorter time interval. From p.62p and (jl.66p we see that 

(1.68) C{t,xo) =Cu,gmt <0, 

so corner compatibility conditions are satisfied in (jl.67p to infinite order. 

We now apply Theorem 11.251 to obtain approximate and exact solutions and to 
the viscous problem (jl.6ip satisfying the estimates (jl.57p for constants eo) M, k, and sq as 
in that Theorem. Note that ti^, u^, and are all equal to u in t < 0. Finally notice that 
smooth solutions ti^ to the initial boundary value problem (|1.6ip are uniquely determined 
from the start. They must therefore equal the solutions obtained from Theorem 1 1 . 25 1 bv the 
above procedure. The estimates involving in ()1.57p now allow us to identify the unique 

solution to (jl.67p as the small viscosity limit of the u^. 
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1.8 Discussion and open problems 

Theorem I1.3UI generalizes to the case of "real" , or partially parabolic viscosity, the small- 
amplitude results obtained in |GG| for Laplacian second-order terms, and to multi-dimensions 
those obtained in [R3j for one-dimensional symmetric-dissipative systems. It includes as a 
physical application existence of small-amplitude boundary-layer solutions for the equations 
of compressible gas dynamics with specified in- or outflow velocity, temperature, and, in the 
inflow case, specified density or pressure, a result analogous to those obtained in |TW| for 
the incompressible case. It includes also the corresponding result for the compressible MHD 
equations with specified inflow or outflow velocity, temperature, magnetic field, and, in the 
inflow case, density or pressure, for parameter regimes satisfying the structural conditions 
described in Remark 15.61 See Sections 15. H 15.21 and 15.31 for physical examples. For a general 
physical survey of boundary-layer behavior, see [Sj. 

Together with pointwise Green function analyses |YZl INZj showing that uniform Evans 
stability is sufficient for long-time (i.e., time-asymptotic) stability of planar boundary-layer 
profiles. Theorem 11.281 yields also long-time stability of small-amplitude planar profiles 
with sharp pointwise rates of decay, sharpening previous results obtained by Rousset |R3j 
by energy methods. 

Theorem I 1 .251 generalizes to the case of real viscosity the large-amplitude results obtained 
in [MZlj and [GRj for strictly parabolic viscosities in the multi- and one-dimensional case, 
respectively, giving a sharp criterion for existence and stability of boundary-layer solutions 
in the small- viscosity limit. Determination of Evans stability in the large-amplitude case 
is an outstanding open problem. (The uniform Evans condition may fail in general for 
large-amplitude layers, as demonstrated in |SZ| .) Numerical testing of the Evans condition 
for large-amplitude layers in multi-dimensions would be an interesting direction for further 
investigation; see [CHNZl HLyZlj for the one-dimensional case, [HLyZ2] for the multi- 



dimensional shock case. 

We note that small-amplitude stability for symmetric-dissipative systems might be prov- 
able for variable-multiplicity systems under weaker structural assumptions than those of 
Theorem 11.301 which are tailored for large-amplitude layers, by direct energy estimates as 
in [GGl IR3] rather than by first passing to the constant-layer limit. This approach becomes 
quite complicated in the multidimensional case, but would yield existence without hard- 
to- verify structural conditions, and would apply to some physical cases such as MHD for 
parameter regimes different from that described in Remark 15.61 

Finally, we discuss the meaning of the somewhat unexpected instability result of Propo- 
sitions [LT3] and Corollary 11.291 parts (c) and (d). In cases where N" > N'^, or equiv- 
alently, when the number of scalar Dirichlet conditions imposed in the viscous problem 
{N' — N" -\- N^) is strictly less than the number of scalar boundary conditions for the resid- 
ual hyperbolic problem (A^+), this seems to indicate a failure of our basic Ansatz, which 
assumes that the residual hyperbolic problems should involve only Dirichlet boundary con- 
ditions. For example, in the extreme case N = N' (strict parabolicity) , = N = Nf^ 
(so A'^^ = and all characteristics are incoming), with full Neumann boundary conditions 
duu(0) = on dQ (so A^" = A^), work in progress indicates that solutions converge in the 
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small-viscosity limit to solutions of the inviscid problem Co{u^) = with Neumann bound- 
ary conditions di^u{0) = on , and with no intervening boundary layer. In other cases 
where A^" > N'^, we conjecture that the correct model for limiting behavior is a residual 
hyperbolic problem with mixed Dirichlet-Neumann conditions of appropriate ranks. We 
plan to address this issue in a future work. 

Plan of the paper. In Section [2] we investigate the existence and transversality of 
boundary layers, in particular in the small-amplitude limit, and construct both small and 
large amplitude C-manifolds. In Section[3l we review the construction of the Evans function, 
and examine its high-frequency and small-amplitude limits, establishing the key reduction 
of Theorem 11.281 for small-amplitude profiles. In Section HJ we establish Corollary 11.291 for 
symmetric dissipative systems using energy estimates. In Section 14.31 digress slightly 
to show maximal dissipativity of hyperbolic boundary conditions associated with small- 
amplitude layers of symmetric dissipative systems when full Dirichlet conditions are imposed 
in the viscous problem (N" = 0). In Sections I5.1H5.31 we carry out explicit computations 
for the example systems of isentropic gas dynamics, full gas dynamics, and MHD. The 
construction of approximate solutions is presented in Appendix [Aj The tracking lemma 
and its connection to construction of high-frequency symmetrizers are given in Appendix 

m 

Notations 1.31. We do not distinguish between and U2, J* and J^, etc.. Sometimes, 
especially when other subscripts or superscripts are involved, one choice is more convenient 
than the other (e.g., u1). 

2 Existence of C-manifolds and layer profiles 

In this section we will show that it is possible to construct smooth C-manifolds as in As- 
sumption [TTT21 globally defined on d^l with corresponding global smooth families of small- 
amplitude layer profiles. We will also show that it is possible to give a local verification 
of Assumption 11.121 if one starts with a given, possibly large amplitude, transversal layer 
profile. 

2.1 Global C-manifolds and families of profiles in the small-amplitude 

case 

Here we verify Assumption 11.121 in the small-amplitude case for a variety of boundary 
conditions of the form ()1.4p . We begin by defining a family of constant layers by giving a 
smooth presciption of states p{t, xq) G with 

(2.1) {t,xo,p{t,xo)) e Uq for all {t,XQ) G [-T,T] x dVL. 

We shall not try to show that such smooth global prescriptions are always possible under 
structural Assumption ll.il but they clearly exist for the Navier-Stokes and MHD equations 
(see sections 15.11 15.21 and 15. 3p , where one has a simple characterization of the domains 
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of hyperbolicity and non-characteristicity in terms of physical quantities like pressure and 
velocity. The prescription (j2.ip determines a choice of boundary data in (jl.4p or (jl.l6p . 
namely 

(2.2) {gi{t,xo),g2{t,xo),0) := {Ti{p\t, xo)),r2{p\t, xo)),0). 

Remark 2.1. When is a half-space, so that v is constant, prescription (j2.ip is trivially 
constructed, consisting of a single state p. 

Given p{t,xo) as in (j2.ip . define the compact set 

(2.3) B := {{iy{xo),p{t,xo)) : {t,xo) G [-T,T] x dfl} C S'^-^ x irUa- 

Fix (z£, p) € B. The next Proposition, which is a slight modification of Proposition 5.3.5 in 
|Met4j . characterizes all possible "small-amplitude" solutions w{z) of 

A,y{w)dzW - dz{By{w)dzw) = on z > 0, 
w[z) q as z ^ oo 

for {i^jq) near (i^, p). Define the N' x A^' matrix 

(2.5) G.{q) := (i?^^)-! (^^2 _ AfiAl^^Al') (q) 

and let Eip(G;y(g')) denote the generalized eigenspace of G^iq) associated to eigenvalues ^ 
with ib3f?/i < 0. Denote by Ii.,y±{q) the projections associated to the decomposition 

(2.6) M^' = E+(G,(g)) e E^{G,{q)), 

and fix isomorphisms a(z/, g; a) linear in a G E_(Gi/(p)) and in (j^, g): 

(2.7) a(i/,g;a) :E_(G,(^)) ^E_(G,(g)) 
such that Q!(z^,p; a) = a. 

Proposition 2.2. There exists a neighborhood oj C S'^^^ x o/ (i£,p) and constants 
R > 0, r > such that for (z^, q) E lo, all solutions w of (12. 4p satisfying 

(2.8) ||(ti;,z/;2)_(^^0)||^^^^^^^ 

are parametrized by a C°° function w = $(z,z/, g, a) on [0, cxd) x iv* , where uj* is the set of 
{v,q,a) with {y^q) G u; and a E E_(G,/(p)) wzi/i |a| < r. T/ie function ^{z,v,q,a) is the 
unique solution of (12. 4p satisfying the boundary condition 



(2.9) U^^w'iiO) = a{u,q;a), 
and /ias i/ie expansion 

(2.10) I/, g, a)=q^ + e'^'^^'>^G-\q)a{u, q; a) + 0(|ap) 
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uniformly with respect to [z,q,v). Moreover, there exist positive constants 5 and C such 
that for all z G [0, oo) and (q, a) G uo* : 

(2.11) \d,^^{z, V, g, a)| + |$(z, V, q, a) - q\ < Ce-^\ 

We also denote by $(z,z^, g, a) the maximal extension of ^ to z < as a solution of (j2.4p . 

Proof. 1. First we claim there exists a neighborhood lo' C S'^~^ x of (i£, p) and constants 
i?' > 0, r' > such that for {v, q) £ uj' , all solutions w of (j2.4p satisfying 

(2.12) lk^llLi[o,oo] < R' and ||w^||^^[o^^] < R' , 

are parametrized by a C°° function w = <I>(z, v, q, a) on [0, oo) x lv'^,, where lo'^, is the set of 
{u,q,a) with {v^q) G and a G IE_(Gy(p)) with |a| < r' . This may be established by a 
contraction mapping argument identical to that given in Proposition 2.2 and Appendix A 
(both) of [GMWZ7] . This argument corrects a minor error in |Met4] . Proposition 5.3.5 and 
extends that Proposition to the case of partial viscosity. 

2. For V near z^, after shrinking and r' if necessary and renaming as and r, the 
map (g, a) (<I>, i^, g, a) defined for (i/, q) G a;, |a| < r, defines a diffeomorphism onto 

the local center-stable manifold of (p, 0) for (j2.4p considered as a first-order system (jl.lSp . 
On the other hand for v near v_ all solutions of (j2.4p for which 

(2.13) \\{w,wl) - (p,0)||ioo[o,oo) is small 

lie on that center-stable manifold. Thus, for R small enough the assertion in the Proposition 
holds. □ 

Remark 2.3. Alternatively, the result of Proposition \2.2\ may be obtained directly by in- 
variant manifold theory, working with the first-order system (jl.lSp . For, noting that any 
constant function is an equilibrium of the system, and recalling that equilibria lie on any cen- 
ter manifold, we find by a dimensional count that the center manifold of the system consists 
entirely of equilibria, and the center-stable manifold is foliated by the union of stable man- 
ifolds through each equilibrium (constant state). Thus, the only profiles satisfying (j2.8p are 
those lying on stable manifolds of rest points q, whence we obtain both the parametrization 
by ^ and the decay estimate (|2.1ip by an application of the Stable Manifold Theorem. 

The next Proposition gives necessary and sufficient conditions on the boundary condi- 
tions T{w,0,dzw'^) in (jl.l6p for the local (with respect to (i^, p)) existence of transversal 
profiles and corresponding C manifolds. The local objects will then be patched together 
using an argument based on local uniqueness. 

For {v,q,a,p) near (i^, p, 0,p), define 

(2.14) *(z., q, a,p) := T($(0, u, q, a), 0, d,^\0, v, q, a)) - (Ti(/), T2(p2), 0) 
Observe that ^'(z^,p, 0,p) = and that every solution of (j2.4p which also satisfies 

(2.15) T(«;,0,w;2)(o) = (Ti(/), T2(/), 0) 
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corresponds to a solution of ^{i', q, a,p) = and vice versa. Using (jl.lSp and the expansion 
(j2.10p . we readily compute the Nb x {N + N') derivative matrix 



/Tl(pi) T[{p') {-{Al')-\p)A'^{p)G^Hp))\ 
(2.16) V,,2,Ji.,p,0,p)= T,ip^)G^Hp) 

V I 

where, for example, the matrix entries in the third column, reading down, have sizes x 
N', {N' - N") X N', and N" x N' respectively. 

Proposition 2.4. (a) For B as in ()2.3p . let {v_^p) G B. The constant layer ^{z,u,p,0) =p 
is transversal if and only if 

(i)the Nh X N' third column of (j2.16p is injective on E,-{Gy{p)), and 
^'^'^^'^ (ii) ifN" > 0, is of full rank N" on E_(G^(p)). 

(b) Suppose ()2.17p holds. There is a neighborhood lu C S'^^^ x -kUq of {v_,p) and for each 
E uj, there is a manifold Ci,,p of dimension N — and a smooth map 

(2.18) Wi,^p : [0, oo) X C^^p ttUq, 

such that for each q G Cu,p, Wy^p{-,q) satisfies (j2.4p . (|2.15p and converges at an exponential 
rate to q as z ^ oo. Moreover, the manifolds Cp^p vary smoothly with {y,p) G a;. 

(c)The endstate-manifolds Cy^p and profiles Wy^p{-,q) are uniquely determined by this 
construction for {q^v^p) near {p,]y,p). More precisely, when (z/, p) lies in charts centered 
at two distinct base points (Ek^P/J ^ B, k = 1,2, the corresponding manifolds C^p are the 
same near p G C^p n C^p and for each q G C^p n C^p, the profiles w^p{-,q) constructed in 
the separate charts coincide. 

Proof, (a). The first transversality condition in Definition 11.111 is equivalent to injectivity 
of ^'a(l£,P, 0,p) on E_(G;/(p)), while the second transversality condition there is equivalent 
to surjectivity of 

(2.19) ^',,a(i£,P,0,p) : X E_(G,(p)) ^ M^" 

(for more detail see |Met4j . Prop. 5.5.3). Condition (i) in ()2.17p is equivalent to the first of 
these conditions, and in view of Assumption 11.71 condition (ii) is equivalent to the second. 

(b). Since ^'0(1^5^1 0,p) has rank N"^ on E_(G,y(p)) and '^q^a as in (j2.19p has rank 
Nb = Nj^ + A'^^ , the Implicit Function Theorem implies that there exist smooth functions 
q{q~, i^,p), a{q-, y,p), where g_ G R^~^+ is a vector consisting of N — N+ of the coordinates 
of q, such that the solutions of ^{v,q,a,p) = near (i^, p, 0,p) are given precisely by 

(2.20) {u,q{q-,v,p),a{q-,u,p),p) for {q-,v,p) near {p_,u,p). 



26 



For each (v^p) the manifold C^^p is defined hy q = q{q-, f^p) and the profiles are given by 

(2.21) Wu,p{z,q) = ^{z,iy,q{q_,iy,p),a{q^,u,p)). 

(c). Suppose {vjP) lies in charts centered at two different base points (EkiPf^) ^ B, 
k = 1,2. Let C^p and w^p denote the corresponding manifolds and profiles. The properties 
of the functions $ and ^ described in Proposition 12.21 and the discussion following 12.141 
show that each of C^p, k = 1,2 coincides near p with the set of q such that there exists a 
w{z) satisfying ([23]),' (f2J5|) . and 

(2.22) \\{w,wl) - (g, 0)11^00(0,00] is smah . 

This description is chart-independent so the manifolds must agree near p. 

Suppose q G C^ p Pi C^p. Working in chart 2 and using the properties of the functions <I> 
and ^ just referred to, we conclude that any small amplitude profile satisfying (j2.4p . (j2.15p . 
and (12.22P must be given by wlp{z, q). In particular, we must have wlp{z, q) = w1p{z, q). 

□ 

Corollary 2.5. 1. For pure Dirichlet conditions (N" = rankTs = 0), sufficiently small- 
amplitude layers are transversal, and likewise for mixed Dirichlet-Neumann conditions in 
the extreme case rankTa = N"^ = N' . 

2. When the number of Neumann boundary conditions rankTs exceeds N"^, constant 
layers are non-transversal. Equivalently, a necessary condition for transversality of constant 
layers is that the number of (scalar) Dirichlet conditions for the parabolic problem, {N' — 
N") + A^^, is greater than or equal to the number of Dirichlet conditions for the residual 
hyperbolic problem, N^. 

Proof. When N" = 0, the condition (j2.17p (i) holds since T2{p'^)G~^ (p) is an invertible 
N' X A^' matrix. When N" = N' = N'^, Assumption [TTl together with the fact that 
dimE_(Giy(p)) = A'"^ imply that is invertible on K^{G,^{p)). Thus, both conditions in 
(pTTl) hold." When A^" > A^^^ it is impossible for plT]) (iilto hold. The final assertion 
follows by noting 

Nh = N" + (A^' - A^") + nI = N+ + N'^ so 
(2 23) 7 I + + I - 

A^" < A^^ ^ {N' - N") +Nl> N+. 

□ 

Proposition 2.6. Given a smooth global assignment of states 

(2.24) {t,xo,p{t,xo)) G Uq for all {t,xo) e [-T,T] x dQ, 

suppose that transversality (equivalently, condition ()2.17p of Proposition \2^)^ holds at every 
point of 

(2.25) B := {{u{xo),p{t,xo)) : {t,xo) G [-T,T] x dft} C S'^-^ x nUa. 
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Then for boundary data 

(2.26) {gi{t,xo),g2{t,xo),0) := {ri{p'{t,xo)),r2{p\t,xo)),0) , 

there exists a global C manifold and family of profiles W satisfying the conditions of As- 
sumption \1.1SX 

Proof. Use the compactness of B to cover B with a finite number of open patches centered 
at basepoints p)a;- Carry out the construction of Proposition YTM part (b), in each patch 
and use the local uniqueness described in part (c) to define manifolds C^^p and profiles 
w^^p{-,q) for all (i^, p) in the covering and smoothly varying with {v^p). Then the conditions 
of Assumption 11.121 are satisfied by taking 

27) Cu(xo),p(t,xo) 

W{z,t,Xo,q) := ^yi.(xo),p(t,xo)(^;9)- 

□ 

Remark 2.7. 1. Part (a) of Proposition can be used to produce many examples of 
transversal constant layers involving mixed Dirichlet-Neumann conditions for any triple 
{N",Nl,N') satisfying < N" < N"^ < N' . Together with Proposition this yields 
examples where Assumption is satisfied for the Navier-Stokes and MHD systems with 
a variety of boundary conditions. 

2. In the case rankTa = N"^ = N' , the matrix Ky is an invertible N' x A^' matrix, so 
the T3 boundary condition in (I1.16P is equivalent to ^^(0) = 0. With (ll.lSp this shows that 
the only solutions of the profile ODE ()1.16p are constant layers. Since N' — N" = 0, T2 is 
absent and we have 

(2.28) C{t,XQ) = {q:Ti{q)=g\t,XQ)}. 

Thus, the residual hyperbolic boundary conditions in this case are the same as the Dirichlet 
boundary conditions for the parabolic problem, with Neumann conditions ignored. The lead- 
ing term in the approximate solution (11.350 is given bylA^{t,x, |) = u^{t,x), the hyperbolic 
solution with no intervening boundary layer. 

2.2 Local C-manifold associated to a given transversal profile 

Proposition 2.8. Suppose that w is a given, not necessarily small amplitude, transversal 
layer profile satisfying 

(a) Ay{w)dzW — dz{BiAw)dzw) = on z > 

(2.29) (6) T{w,0,wl){0) = {g\g^O) := {g,0) 

(c) w{z) p as z ^ 00. 
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Then in a neighborhood O C of p and for {v,g) near {]y,g), there is a smooth manifold 
Cu,g C O of dimension N — and a smooth map 

(2.30) W^^g : [0,oo) xC^^g^TTU*, 

such that for each q E C^^g, Wu^g{-,q) satisfies l\2.29^ (a).(b) with {i^jg^q) in place of {u,g,p) 
and converges at an exponential rate to q as z oo. Moreover, these manifolds and profiles 
vary smoothly as {v^g) varies near {]y,g). 

Proof. 1. For $(z, z^, q, a) and 0(1^, q; a) as in Proposition 12.21 the given profile must satisfy 
for some sufficiently large zq 

(2.31) ^(2:) = ^{z — zo,i/,p,a) for a such that U^^^vuI^zq) = a{iy,p;a). 

This follows from the fact that for some zq sufficiently large, the condition (j2.8p of that 
Proposition holds with [0, 00) replaced by [zq, 00) (see Prop. 5.3.6 of |Met4] for details). 
2. For {i',q.,a,g) near (i£,p, a, 5), instead of (j2.14p we now define 

(2.32) ^{v, q, a, g) := T^O, q, a), 0, d,<^\0, u, q, a)) - {g^ , 5', 0) 

and observe that ^{E.^p,a-,g) = 0. Transversality of w implies that 

^a{R,P,a,g) : E_(G,(p)) ^ M^" and 
*g,a(id,P,a,9) : x ^~{GJj))) ^ M^'' 

have ranks A''^ and = + N'^ respectively. The Implicit Function theorem implies that 
there exist smooth, locally unique, functions q{q-,i',g), a{q-,u,g), where g_ S M^~^+ is a 
vector consisting of N—N^ of the coordinates of q, such that the solutions of ^{i', q, a,g) =0 
near (i^, p, a, g) are given precisely by 

(2.34) {u,q{q^,u,g),a{q^,iy,g),g) for {q-,v,g) near {p-,E,g)- 

For each (z^, 5) the manifold Ci^^g is defined hy q = q{q-, u,g) and the profiles are given by 

(2.35) W^,g{z,q) = ^>(z - zo,u,q{q-,iy,g),a{q-,iy,g)). 

□ 

2.3 Global C-manifold for a family of transversal profiles 



Similarly as in (j2.ip for the small-amplitude case, assume that we are given a smooth 
prescription of large-amplitude profiles, in the form of a C°° function 



(2.36) w{z, t, xq) : [0,00) X [-T,T]x dn 
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that defines a transversal layer profile for each (t,xo): 

(a) A^^^^^){w)d^w - d^{B^^^^^-jdzw) = 
^2 (&) 0' wl)\z=o = {^\wH0, t, xo)), T\w\0, t, xq)), 0) 

:= (5^(^,2:0), 5^(i,a;o),0) 
(c) ^ ty(oo,t,j;o) := q{t,xo) as z ^ 00. 

Remark 2.9. When il. is a half-space, so that v is constant, assignment (j2.24p exists 
trivially, consisting of a single profile. 

Corollary 2.10. Given a smooth transversal family ()2.36p . there is a smooth manifold C 
defined as the graph 

(2.38) C = {{t, xo,C{t, xo)) : {t, xq) e [-T, T] x OQ} C U9, 



where C{t,Xo) C M is an N — dimensional manifold containing q{t,xo) and consisting 
of states r near q{t,XQ) for which there exists a transversal layer profile 

(2.39) W{z,t,XQ,r) : [0,00) X C ^ nW 
satisfying (12.37^ with g(t,xo) in (c) replaced by r and 

(2.40) \\W{z,t,xo,r) - w{z,t,xo),dz{W'^{z,t,xo,r) - uj'^{z,t,xo))\\L^ small . 

The profiles W{z,t,X(),r) are C°° in all arguments and converge to their endstates at 
an exponential rate that can be taken uniform on compact subsets ofC. 

Proof. 1. Cover the compact set B = {{v^xq), q{t, xq)) : (t,xo) G [— T, T] x dil} by a finite 
number of charts centered at points (z^, q)k for which we have defined 

(2.41) a{iy,q;a):E-{Ciqk))^E-{G,{q)) 

as in ()2.7p . Fix (t,xo). As in the proof of Proposition 2.8, we have for some k 

(2.42) w{z,t,xo) = <^>{z-T,u{xo),q{t,Xo),a), a £ {G^^iqk)) 
for some large T, where a is such that 

(2.43) n^(a,g) _(g(t,xo)) dzw'^{T,t,xo) = a{v{xo),q{t,xo),a) e E~{G^(^^^j{q{t,xo))), 

with Il^^^{q) the projection of onto {Gu{q)) along E+(Gjy(g)). 

2. Setting W{z) = {w{z),w'^{z)) we rewrite (|2.37p fa) as a first-order system: 

(2.44) d,w = n,^,,-jiw). 
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rr ^\\'-^^^ near a} 



The stable manifold of xo),0) for ()2.44p near {w{z,t,xo),w'l{z,t,XQ)) for any fixed 
z G [0, oo), can be parametrized: 

w%q{t,xoy,iy{xoy,z) = 

(2.45) J^z-T,u{xo),q{t,Xo),a) 
^<^l{z - T,i'{xo),q{t,xo),a) 

Similarly, the center-stable manifold of (g(t, xo),0) near {w(z,t,xo),'Wz{z,t,xo)) is 

W'''{q{t,xo);i^{xoy,z) = 

^^•^^^ ^ (<^t - ^r^^tlt a)) ^ ^ ^'(GuM) near a}. 

If z = T we get the parts of the respective manifolds near 

{w{T,t,xo),w'i{T,t,xo)). 

If z = we get the parts near {w{0,t, XQ),wl{0,t, xq)). (Here we use that $(z,z^, g, a) is 
a maximal extension to z < 0). Note that for any fixed z, these manifolds are uniquely 
determined locally near {w{z,t, xo),w'^{z,t, xq)). 

3. The set 

(2.47) {(r,0) :r GC(t,xo)} 

is, near {w{oo, t,xo),0), the set of endstates at infinity under the flow of (j2.44p of the N—N^ 
dimensional initial manifold 

(2.48) {U = {u\u\u^) : {Ti{u^),T2{u^),0) = {g\t,xo),g\t,xo),0)}n 

>V^^((z(t,xo);K^o);0)cM^+^', 

where we have evaluated (I2.46P at z = 0. As a check observe that the intersection (12.48P is 
transversal and has dimension 

(2.49) {N + N' - Nb) + iN + nI) - {N + N') = N - N+. 
Considering the uniqueness of W^^{q{t, xq); z^(xo); 0) near 

{w{0,t,xo),w1{0,t,xo)), 

this description of C(t, xq) establishes that it is uniquely determined by the local construction 
of Proposition 2.8 near w{oo,t, xq). This allows the global C manifold of Corollary 12.101 to 
be constructed by patching together local C manifolds. 

4. The profiles W{z,t,xo,r) in ()2.39p are given by 

(2.50) W{z, t, xo,r) = Wi,(^^^)^g(^t,xo){z, r) 

for Wp^g as in (2.33) of GMWZ5. For each point (r, 0) in (j2.47p . there is a unique point Uq = 
{uq,Uq) G Cinitiai{t,XQ) that is mapped to it under the flow of ()2.44p . w{z) = W{z,t, xo,r) 
is uniquely characterized by the properties that: 
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a){w,Wz) satisfies (|2.44p 

c) {w,w'^) — i- (r, 0) as z — > +00. 

d) {w,w'^J{0) = Uo £ Cinitial{t,Xo). 

This characterization shows that the locahy constructed profiles 

W{z,t,xo,r) 

in (j2.50p patch together consistently. The local construction of Prop. 2.8 then shows that 
they are C°° in all arguments. 

□ 

3 The Evans function and asymptotic limits 

We now study the Evans function and its high-frequency and small-amplitude limits. We 
first recall the conjugation lemma of |MZ1] . which implies that a first-order system U' = 
Q{zX,p)U whose coefficient matrix converges exponentially to its limit Q{oo,C,p) as z ^ 
+00, may be converted by a smooth, exponentially trivial local change of coordinates 

(3.1) U = P{z,C,p)V = {I + Q{z,C,p))V 

to its limiting constant-coefficient equation V = G{oo, (^,p)V. Here, we have adjoined to the 
arguments z also dependence on model parameters p, assumed to be at least continuous. 

Let G{z,CjP) be as in (|1.45p . a frequency-dependent matrix arising from linearization 
around a profile w{z) such that for some positive constants C, /?, uniform with respect to 
model parameters p, 

(3.2) \w{z) -w{oo)\ <Ce-'^^, 

and also p {w, dzW2){-,p) is continuous as a function from p to L°°{0, +00). Thus, also, 

(3.3) \g{zr)-G{oo,-)\<Ce-''' 
and g{-,C,p) is continuous as a function from p to L°°(0, +00). 

Lemma 3.1 ( [MZT] . Lemma 2.6). Let P > be as in ([3:2]) . For all C G M''+^ and model 
parameters p, there are a neighborhood lo x P of {C,p), « matrix P{z, (,p) = I + Q{z, C,p) 
that is on [0, 00) x uj with derivatives uniformly continuous in p, and positive constants 
C , a with < a < /? such that 



(i) P and P ^ are C°° and bounded with bounded derivatives: 

(3.4) \d{dlQ\<C,ke-''\ 
(a) P{zX,p) satisfies 

(3.5) dzP = g{z,C,p)P-Pg{^,C,p) onz>0. 
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Observe that U satisfies (jl.45p on z > if and only if V defined hy U = P{z,C,p)V 
satisfies 

(3.6) d,v = g{^,c,p)v + p-^F, rp{o,c,p)vu=o = G. 

Tliis implies that the decaying space E^((",p) as in (jl.48p is exactly the image under 
P{0,CiP) of the stable subspace of ^(oo,C,p), denoted K'^{(,p). Thus, by the calcula- 
tion of |GMWZ6j . Lemma 2.12, E-(C,p) has dimension Nf, = rank T for 7 > 0, C / 0. The 
Evans determinant p.48|) 

(3.7) DpiO = det(E-(C,p),kerr(C,p)), 

now denoted with additional dependence on model parameters p, is then well-defined on 
7 > 0, C / and depends smoothly on ( and continuously (in all C-derivatives) on p. We 
record this as a corollary, of which we shall later make important use. For quantitative 
bounds estimating the modulus of continuity in p, see |PZj Prop. 2.4 or Cor. C.3, |HLZ] . 

Remark 3.2. The conjugator P{z,(,p) is constructed by a fixed point argument as the 
solution of an integral equation. The exponential decay ()3.3p is needed to make the in- 
tegral equation contractive in L°°[M, +00) for M sufficiently large. The continuity of P 
with respect to p is then immediate, by continuous dependence on parameters of fixed point 
solutions, a quite general result. 

Corollary 3.3. Let w{z,p) be a family of layer profiles depending continuously on parame- 
ters p in the sense that p ^ {w{-,p),w^^{-,p)) is a continuous function from p to L°°[0, +00), 
and letV{C,,p) be as in (13. 7p . Then, the Evans function (13. 7p depends continuously on p. 

Proof. 1. Set W{z,p) := {w{z,p),uj'^{z,p)). Continuity in p, by boundedness of Aj, Bjk, 
and derivatives, is inherited by the coefficient matrices Q{-, C,p) appearing in the linearized 
eigenvalue equations from continuity of W . Likewise, continuity of the linearized boundary 
operator T{C,,p) follows from boundedness of T and derivatives. In view of our rank condi- 
tions on T{C,,p) and the continuity of -P(0, C,p) for all C, G M'^"''^, we see from the definition of 
the Evans function (j3.7p that it is sufficient to establish continuity of Ej^((^,p) with respect 
to p for C / and continuity of E;^(C, /9,p) with respect to p at /? = (recall (jl.5ip . (|1.52p ). 

2. For C 7^ the continuity of E^(^,p) follows by the fact that the limiting coefficient 
matrix Goo{C-,p) has a spectral gap, whence the stable subspace varies continuously by 
standard matrix perturbation theory |Kat j . Continuity Ej^(^, p,p) at p = is more difficult 
to show and follows from the existence of ii'- families of viscous symmetrizers as in ^MZ3] . 
Theorem 3.3. 

□ 

3.1 Maximal stability estimates and high-frequency scaling 

We next recall from |GMWZ6] the appropriate scaling of the Evans function for high- 
frequencies I CI > i? > 0. The maximal stability estimate for (jl.45p on this frequency 
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domain is (see |GMWZ6] ) 

(1 + 7)II'"^IIl2{m+)+A||u^||l2(]r^) + ||9^n2||L2(R+) 
^3g^ +(l + 7)^|ni(0)|+A5|n2(0)|+A-5|5.n2(0)| < 

C(II/^IIl2(m+) + A"-^||/^||l2(]r^)) 
+ C{{l+j)-2\g'\+A"^\g^\+A'-2\g^\), 

where C is an independent constant and A is the natural parabohc weight 

(3.9) MC) = {r'+l' + \v\Y'- 

Together with corresponding low- frequency estimates (see (I3.16|] below), this yields maximal 
spatio-temporal stability estimates by Parseval's identity [GMWZ4[ IGMWZ6| . 
Taking / = in l\3.8\i yields the necessary condition 

(3.10) (1 + 7)^|m^| + A^|u2)| A-3|u2| < 

C((l + j)"^\Tiu'\ + A-2\T2U^)\ + A--2mO{u\ u^)\) 

VC G 1+ , Id > R,yU = (n\u2,u3) G E_(C). This can be reformulated in terms of a 
resettled Evttns funetion (see [MZl] ). Introduce maps defined on C^"*"^' and C^*" respectively 
by 

^^^^^ Jdu\u\u') := {{1+j)'2u\a'2u\A--2u') 

Note that J(r{C)U = r'^{C)J^U with 

(3.12) = (TmS ^2M^ KdU^ + a-^Kt{v)u^) 
(note: decoupled, bounded). Thus (j3.10p reads 

(3.13) VC/ E JcIE-(C) : \U\ < C\J^T{C)J-^U\ = ClT'^Ul. 
Introducing the resettled Evttns funetion 

(3.14) D^'^iC) := |det(jfE-(C),Jckerr(C))| = | det (j^E-(C),ker r^(C)) |, 

and using Lemma II. 2H we see that this stability condition is equivalent to the following 
definition. 
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Definition 3.4. (a) Given a profile w, the linearized equation (jl.42p satisfies the uniform 
Evans condition for high frequencies when there are c > and R > such that \D^^(Q\ > c 

for all C G ij*'^ with |C| > R. 

(b) The linearized equation (jl.42p satisfies the uniform Evans condition when there are 
c > and R > such that 



For completeness, we recah also the maximal stability estimates for low- and medium- 
frequencies Id < R, of 



where = (7 + |Cp)2, for ( £ M+ \{0}, |C| < R- Taking / = yields the necessary 
condition |ti(0)| + |9^u^(0)| < C\g\ corresponding to the standard Evans condition (jl.48p . 

Remark 3.5. By Theorem 3.9 of IGMWZG^ the uniform Evans condition for low and 
medium frequencies implies the maximal estimate (13.16^ . Taking f = in (I3.16P and using 
Remark \1.2S^ we see that the following are equivalent for |C| < R: 

a) the uniform Evans condition for |C| < R, 

h) the full hounded frequency estimate ()3.16p . 

c) the trace estimate in (j3.16p when / = 0. 

Similarly, by Theorem 7.2 of 'GMWZ6^ the uniform Evans condition for high frequencies 
implies the maximal estimate p.Sp . Using Remark \1.22\ again, we deduce the equivalence 



a) the uniform Evans condition for \(^\ > R, 

b) the full high frequency estimate (I3.8p . 

c) the trace estimate in (j3.8p when / = 0. 

3.2 The high-frequency hmit 

In this section we show that the rather complicated high-frequency condition of Definition 
13.41 mav be reduced to a simple and natural linear algebraic condition corresponding roughly 
to well-posedness of the principal part ()1.8|) of equations ()l.ip . frozen at ^; = 0, under 
boundary conditions (jl.4p . Using the fact that the linearized boundary conditions are fully 
decoupled, we prove in Theorem 13.61 that satisfaction of the uniform Evans condition for 
sufficiently high frequencies is equivalent to satisfaction of the uniform Evans condition for 
the decoupled system 



(3.15) 



D{C)\ > c for Id < R and \D'^{C)\ > c for \C\ > R. 



(3.16) 



v\MlHr+) + \\dzU^\\L^(R+) + \u{0)\ + \dzu^{0)\ < 

C{-\\f\\LHR^) + \9\) 



of: 



(a) ni + J]A"(u;(0))a,^xi = 0, 



(3.17) 



j 



(b) ui-Y,B%[wmd,dku^ = ^. 
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with boundary conditions Fi and (r2,r3), respectively. 

Let IE~(C) denote as before the set of initial data at z = of decaying solutions of 
dzU — Q{z,()U = 0. Our proof of Theorem 13.61 is based on showing that in the limit as 
Id — > oo, the space J^E~(C) approaches, or "tracks", the rescaled stable subspace of an 
appropriate frozen coefficient problem. 

Since subsystems (|3.17p (a)-(b) are quasi- homogeneous, the uniform Evans condition is 
equivalent simply to nonvanishing of the decoupled Evans functions on the unit sphere in 
7 > 0, a linear algebraic computation. This reduction to a compact set of frequencies 
(Corollary 13. 7p is essential in our later verifications of high frequency stability. 

To state these results precisely we write the first order systems obtained from (j3.17p (a) , 
(b) by Fourier-Laplace transform as 



{b)d,U* - G2{0U* = 0, T,U* = (5',0), 



(3-18) ^ _ _ ,2 



where U* = {u^^uD and = (F2,F3). Let e^^hiC) aiid ^-,p{C) denote the stable subspaces 
of the matrices Gi{C) aiid ^2(C) respectively. Setting 

Ji(ni) = (l+7)5n\ Ffn^ = Fini 

(3.19) J^*(n2,n3):= (Ain2,A-^n3 
Ff(C)(^x^u3) = {T2U^KdU^ + A-'KT{r])u^) , 

we define rescaled Evans functions for (I3.18p (a) and (b) by 

Df{C)= det (J>e_,,,(C),kerFi) (= det (e_,;,(C), ker Fi) clearly), 

pN — N' ^ 

(3.20) . 

■■= det (j^*e_,p(C),kerFf(C) 

The main result of this section is the following theorem. 

Theorem 3.6. Let D^^ be the rescaled Evans function defined in ()3.14p . Then frozen- 
coefficients stability for (j3.18p . that is, the existence of positive constants ci,C2, R' such 
that 

(3.21) l^f (01 > ci and |I)f (C)| > C2 for all \C\ > R' , 

is equivalent to the rescaled uniform Evans condition, that is, the existence of positive con- 
stants c, R such that 

(3.22) \D'%C)\>c for all\C\>R. 

Proof. 1. Frequency zones. There are two frequency zones to consider. For 5 > 
sufficiently small we define the elliptic zone 

(3.23) ^5 :={C=(t, 7, ?7) :7>5|C|,|r?| >5|C|} 
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and the complementary coupling zone 



(3.24) Cs:={C:0<l< 5|C|} U {C : |r?| < 5|C|}. 

Here we show how the discussion of [GMWZ6] . section 7 can be adapted to prove Theorem 
13.61 in the more difficult case where G C^. The elliptic case is proved in Appendix B, which 
also includes a new discussion (see Proposition IB. 4p of the connections between symmetrizer 
estimates and tracking. 

2. Symbols. Let r"* denote the space of homogeneous symbols of order m, that is, C°° 
functions h{z, Q such that for all a G N'^"''^, all A; € N, and some 9 > 0, there are constants 
Ca^k such that for \(\ > 1, 



(3.25) 



<C«,o|Cr-'°', if ^ = 0, 
\d^d^h\ < C«,fee-^^|Cr"l"l, if A; > 0. 



Let PF"^ denote the space of parabolic symbols of order m, that is, C°° functions h(z, C) 
satisfying similar estimates with \(\ replaced by A. 

3. Conjugation to block diagonal form in Cg- Consider again the linearized 
problem (jl.45p 

(3.26) d,U = giz,OU + F, T{C)UiO) = G. 

where the components of G are given explicitly in (|B.33p . Set 



(3.27) V 



22 



IC|/ 

|^|-lg32 g33 



In Lemma 7.3 of |GMWZ6] it is shown that there exist positive constants c, R, and 6 such 
the distance between the spectrum of G^^{z, (^) and the spectrum of V'^'^{z, () is larger than 
c|C| for C G Cs, \C\ > R- For such Lemmas 7.5 and 7.6 of [ GMWZ6] use this spectral 
separation to construct a change of variables U = VU, F = VF that transforms (|3.26p to 

(3.28) d,u = g{z,c)u + F, r{Ou{o) = G. 

Here Q = Qp + Q' and with obvious notation 



(3.29) g. 



p 



We have V = Vi{z, ()Vn{z, (), where 

/ 0^ 

(3.30) Vi = I ICr^V^i / 

V^i /> 



gll 













r-1 
















pi 





g32 _ y31gl2 








ro 


po 



yl2 


|C|-ivi3\ 




I 




, V'^ G P° 
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4. Rescaling. It is helpful to work with an equivalent rescaled system. Let J := Jq be 
the scaling operator introduced in (j3.1ip (There are really two such operators, but we use 
J to denote both.) For any matrix M of size {N + N') x {N + N') or Nb x {N + N'), and 
for any U G C^+^' or C/ £ C^" define 

(3.31) Ms = JMJ-\Us = JU. 



In expressions like Us or Qs where the order of the "hat" and scaling operators may be 
unclear, the scaling operator always comes last. 

With V = Us the system (j3.28p is then equivalent to 



(3.32) 



dzV = Gsiz, C)V + V^'Fs, f ,(C)y(0) = Gs. 



Observe that = TsV^ and that Tg is the same as the operator F*"^ defined in (j3.12|) . We 
may write Qs = Ops + G's where 



(3.33) 



'g^^ 

A 

A"l(g32_v31gl2) g33 



F-1 (1 + 7)3A-^F0 (l + 7)iA5F-i' 



A5(l+7)-5F-^ F° 
\A-|(l + 7)-iF0 A-^FO 



AF-i 



It will be important to know the exact form of V^. Direct computation of JVJ gives 



(3.34) 



1 




(1 + 7)2A"2V 
/ 

,-li;31i;12 



12 



\A-2(1 + 7)-2V^i A-^V^^V 



The only entry of this matrix that is not obviously bounded as |C| ^ oo is the (12)-entry. 
By equation (7.36) of [(^IMWZBj we have 

(3.35) V'\z,0 = O{Hm), 

so boundedness of the (12)-entry follows from 

(1 + 7)^ 



(3.36) 



(1 + 7)5|?7|/|CI 



r < |??| 2 < A2 . 



A similar computation of V^^ shows 

(3.37) \Vs\ < C, \Vs^\ < C uniformly for z > 0, C e Cs, 
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and furthermore 

/I 0(1) 0(1)\ /I 0(1) 0(1)\ 

(3.38) V, PS 1 , « 1 for Id large. 

\0 I I \0 1 / 

5. Incoming-outgoing estimates. Consider now the principal part of ()3.32|) : 

(3.39) d,V = gps{z,OV + Vi^Fs, TsVsV{^) = Gs. 

By Lemma 7.13 of [GMWZ6] there is a smooth change of variables W E PT^ such that, if 
we set 



(3.40) r(z,C):= 
then 



/ 
Wj ' 













(3.41) 


T-^QpsT = 


p+ 







V 








with P± having their eigenvalues satisfying ibJR^ > cA. Defining W by 

(3.42) V = TW, W = {w^,w+,w-), 
we can write ()3.39p equivalently as 

(3.43) d,w = gB{z,c)w + T-^v-^F„ r,V,TT^ = o,, 

where now Qb = Gb,p — T~^Tz. 

Define the outgoing (VF^) and incoming {W~) parts of W by 













(3.44) 














K j 







where w^~^ = , =0 when N\ = 0, w^^ = = when N\ = N — N' . With 

ll^ll •= we define norms 

(a)||l^||s = (1 +7)^||tu^|| + A5||t(;+|| + A3||tt;-|| 
(3-45) {h)\\F\\'^ = (1 +7)~5||f1|| + A^^IIF^II + A~3||i?3|| 

{c)\wms = \wm- 

Proposition 7.11 and Corollary 7.14 of |GMWZ6] imply that for large enough R and 
C G with \Q\ > R, we have the following estimates: 

.3^g. \\w+\u + \w+m<c\\{d,-gB,p)w+l 

\\w-\\s<c\\id,-gB,p)w-t + \w~io)\. 
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Suppose now that F is a solution of (j3.32p . Let T be exactly as above and define W by 
V = TW, W = {w'',w+,w-). Now W satisfies 

(3.47) d,W = gBiz,C)W + T-^g',TW + T-^V;^Fs, T,VsTW = G„ 

/,From (|3.46p we deduce the following estimates for W by treating the extra terms T^^Q'^TW 
and T~^TzW as forcing terms: 

.3 \\W+\\s + \W+m < C\\{d, - gB,p)W+t + e{C)\\W\\s 

\\w~\\s<c\\id,-gB,p)w-\\', + \w-io)\+eiO\\w\\s. 

Here we have used the explicit forms of Q'^ (|3.33|) and T. For example, the first row of 
T-'^g',TW contributes (recall dOSil fb)) 

(3.49) {1 + -f)-^\\r-^w^\\ + C\\A-^r\w+ ,w-)\\ + C\\A^T-\w+ ,w-)\\, 

which is a term of the form e(C)||VF||s. 

Defining the outgoing and incoming parts of V by 

(3.50) V+ = V- = TW~, 

we deduce from (I3.48P and (I3.37P the following estimates for solutions to (j3.32p : 

.35^. («) \\V+\\s + \V+m <C\\F,t + e{C)\\V\U 

ib)\\v-\\s<c\\FX + c\v-m+40\\y\\s, 



where e(C) ^ as |C| 00. 

6. The stable subspace of Gps{0,C)- Let Et{C) be the stable subspace of the frozen 
operator GpsiO, Q- ^From ([OT]) we see that when A^j = (resp. Nl = N - N') 

ELiC) = {TiO,C){ Oj :w-€C^'}, 



(3.52) 



(resp., ,ii;!(c) = {r(o,c) I \ -.w'- eC'-^^w- eC"'}), 



^From this and p.44p we conclude that for V as in (j3.50p 

(3.53) V-{0)eE{{C). 

Proposition 7.10 of [(IMWZGj shows that, just like ^i(C) in (I3l8]) fa). the matrix has 
eigenvalues with only positive, respectively negative, real part if a]1 is outgoing {Nl_ = 0), 
respectively incoming (A^^ = N — N'). Moreover, for := as in ()3.19p and Q2 as 



40 



(j3.18p (b). the (32)-entry of Qpg differs from that of J*Q2{0J* ^ by a term that is 0(1). 
These observations imply that for \(^\ large 

(3.54) eUC) is close to j^'^C^'-^''^^^)^ e (^j,^ ° (recall ^M)- 

7. Incoming-outgoing estimates imply tracking. Suppose that U (z) is a decaying 
solution of 

(3.55) d,-giz,c)u = 0. 

Then U{0) G E~(C), or equivalently, 

(3.56) V{0) = JV'^U{0) G JV"^E"(0). 
^From (I3.5ip we easily obtain 

. . \\V-\U<C\V-m+e{C)\\V^s 

\\V+\\s + \V-^iO)\<C\V~m+e{C)\\V+\\s, 

and thus 

(3.58) |y+(0)|<e(C)|y-(0)|. 
Using ([338]) . ([336]) . and (|332]) we obtain 

(3.59) JV-^E-{C) ^ Ei{C) for ( G C^, |C| > R. 
Applying JVJ'^ to ([339]) we find 

(3.60) J1E"(C) ~ JVJ-^eL{C) = VsEL{Q). 

Since has the upper triangular form (|3.38|) as |C| — > oo, we conclude from (|3.6U|) and 
(I33i]) that 

(3.61) JE-(C) ^ (^^'' o'^'^^) © (j*e°p(C)) - ^'^ ^ 

8. Conclusion. Since 

(3.62) kerr, = ^^^^^^^'^ ~ ^ ^ 



/ IkerP 



c I ' 



()3.6ip implies the equivalence in Theorem 13.61 for |^| > R, C ^ Cs- Together with the proof 
for the elliptic zone £s in Appendix B, this concludes the proof of Theorem 13. 6i 

□ 
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As a corollary we obtain the following simple criterion for high-frequencies, 

Corollary 3.7. Under (H1)-(H6), satisfaction of the high-frequency uniform Evans con- 
dition for 1^1 > i?, > sufficiently large, is equivalent to nonvanishing on the positive 
parabolic sphere 7 > 0, |A| + jr/p = 1 (equivalently, the positive unit sphere 7 > 0, |C| = 
of the parabolic Evans function 

(3.63) d\0 := det (e_,p(C), ker Tf (C)) . 
Proof. The system (|3.18p (b) has the form 

(3.64) d.U*=g2m*:=(^l^ 

where M.{C) and ^(C) are quasihomogeneous of degrees two and one respectively. The 
equation (j3.64p can be written equivalently as 

d,{J^U*) = AQiiOJ^U*, where 

(3.65) /O 1- 

^^2(0 ■= [ M A 

This shows that U* G e_^p(C) if and only if J^U* is in the stable subspace of ^2(C)- With 
:= C/A and writing E-{M) for the stable subspace of any matrix M, we thus have 

(3.66) 4e_,p(C) = E.iG^iO) = E^iG^iO) = e-,p(C). 
Since we clearly have 

(3.67) kerrf(C) = kerrf(C), 

it follows from (|3.66p that |D^(C)| > C > for |("| large if and only if (i^(C) is nonvanishing 
on the parabolic unit sphere. 

Since the eigenvalues of were assumed all positive (totally incoming) or all negative 
(totally outgoing), we have either kerFi = {0}, e_^/j(C) = C^~^ or else kerFi = C^~^ , 
e_^/j(C) = {0}. In either case the hyperbolic stability condition is trivially satisfied. □ 

The following result verifies high-frequency stability for many physical cases, including 
the applications we will consider here. 

Proposition 3.8. Consider a layer profile w{z) as in (|1.4ip and the linearized equations 
about w{z) given by (I1.42P (or (11.45p . In the symmetric-dissipative case (Defn. 
the uniform high-frequency Evans condition is satisfied either for full Dirichlet conditions 
rankTa = or full Neumann conditions rankTa = N' on the parabolic variable . 
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Proof. 1. By Corollary 13.71 and Remark 11.221 the uniform high-frequency Evans condition 
in the case of full Dirichlet (resp. full Neumann) boundary conditions on is equivalent 
to the estimate 

(3.68) 1^2(0)1 < C\g\ (resp. \u\0)\ < C\h\) 
for decaying solutions of 

(3.69) k^d j,k^d 

n2(0) = g (resp. li^(O) = h), 

where the constant C in (13.68P is independent of (A, rj) in the positive parabolic unit sphere. 
As in (j3.17p the coefficients in (j3.69p are evaluated at w{0). 

The estimates below are similar, but not identical, to those given in section 14.11 Here 
we highlight the differences and refer to that section for extra detail. We now take 



(3.70) |A| + |r/ 



2 



2. Dirichlet conditions. Taking the real part of the inner product of with (|3.69p . 
we obtain after integration by parts as in (14.70 : 

(3.71) (3f?A + \v\'')\u2\l + \u^\l < C{\g\\u'^m + \v\\g\^). 

Here the last term on the right is a "Carding error" that is explained just below ()4.6p . 
Similarly, taking the real part of the inner product of —U2 with (I3.69p . we obtain as in 

m 

(3.72) i^X + M'Mll + Kli < C ((|A| + |r/|2) \g\\um + |r?||4(0)P). 

The small differences between the estimates here and in section 14.11 reflect the absence of 
the matrices Ak in (|3.69p . 
Using the Sobolev bound 

(3.73) |u^(0)|2 < |n'2|2|4'|2 < Ql^li + 6\u"\l 
we immediately deduce 

(3.74) |n^(0)| < C\g\ 

from (ITtTD . (I372I1 . and (fOHD . 

3. Neumann conditions. Taking inner products as above, but now taking imaginary 
parts in order to estimate |9A||u2|2 and |9A||ti2l25 we obtain after combining estimates: 

(3.75) (|A| + \v\^)\u2\l + I4I2 < CMu'^m + \v\\g\^) 
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and 



(3.76) 



(|A| + \V\')\U2\1 + K\l < C m + |r?|2) \g\\u',m + H\um?)- 



Using the Sobolev inequality 



(3.77) 



and the estimates on U2 and u'2 coming from the first terms on the left in (j3.75p and (j3.76p 
respectively, we obtain 



Remark 3.9. Recalling (see Remark \3.5\) that the Evans condition is equivalent to max- 
imal stability estimates, we find from the reduction to ()3.17p that for decoupled boundary 
conditions the assumed ranks ofTi and (r2,r3) are necessary in order to obtain maximal 
high-frequency stability estimates. For example, specifying density or pressure rather than 
velocity for outgoing isentropic gas-dynamical flow (see Section [5A\) would result in degraded 
stability estimates. 

Note that the above results hold for arbitrary -amplitude layers. 
3.3 The small-amplitude limit 

With Corollary 13 . 71 we are now able to verify the uniform Evans condition for high frequen- 
cies (Definition 13. 4p by reducing to the consideration of a bounded set of frequencies. This 
puts us in position to prove Theorem 11.281 

Proof of Theorem \1.28[ 1. Preliminaries. It is sufficient to show that uniform Evans 
stability of the constant layers w{z, u, v) = u, {u, 1^) € D implies uniform Evans stability for 
sufficiently small amplitude profiles associated to elements of D. (The reverse direction is 
trivial, zero-amplitude being included in the set of small-amplitude profiles.) By compact- 
ness of D, it is sufficient to establish stability of small-amplitude layers in the vicinity of 
the constant layer w{z) = w{z,u,u) = u associated to a single element [u,]/) £ D. Recall 
that e-amplitude profiles w{z,u,v) as in Definition 11.261 satisfv 



for some e > 0. 

2. Parameters. From the assumed Evans stability of w = n, we have transversality of 
the constant layer by Lemma ll.23[ Thus, by Proposition 12.21 specialized to the vicinity of 



(3.78) 



U2m<c\h\. 



□ 



(3.79) 



(a) Au{w)dzW — dz{By{w)dzw) = on z > 

(6) u; ^ n as z ^ 00, 

(c) \\{w,wl) - (u,0)||loo < e, \u-y\<e 
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a single point, there exists a neighborhood uj C x S'^ ^ of {u,u) and constants R > 0, 
r > such that for {i^,q) G uj, all solutions w of (|3.79p (a).(b) satisfying 

- (n,0)||i«>[o,oo] < ^, 

are parametrized by a C°° function w = g, a) on [0,oo) x uj*, where uj* is the set of 

parameters {u, q, a) with (u, q) £ w and a G E_(Giy(p)) with \a\ < r. 

Let Du^uiC) denote the Evans function for (jl.26p corresponding to linearization around 
the constant state u, let Dy^q^aiC) be the Evans function arising from a profile w = ^{z, v, q, a). 
and note that 

(3.80) Du,u = D^,u,o- 

Using similar notation for the Evans function d?{C) defined in ()3.63p we have: 

, . (a)4(C) = det (e_,p,fe(C),kern^fc(C)) , k = {u,v), {v,q,a) 

(6)^fc(C) = det(E^(C),kerrfc(C)), fc = (n,iv), (z.,g,a). 

Observe that for v near v_, given 5 > there exists < e < i? such that 

(3.82) \\{w,w'l) — {u,Q)\\La3yQ^^-^ < e, \u — u\ < e ^ \v — u\ + \q — u\ + \a\ < 5. 

This follows from the fact that by Proposition 12.21 for ^ near u_, (<I>,<I>^) defines a diffeo- 
morphism from a neighborhood of {q, a) = {u, 0) into the center-stable manifold of {u, 0) 
for (j3.79p (a). written as a first-order system. 

2. High frequencies. By Proposition 13.71 the uniform Evans condition for high fre- 
quencies is equivalent to the existence of c > such that |(i^(C)| > c for C, on the positive 
unit sphere 5*+ := {C, : |C| = liT ^ 0}, a compact set. Thus, it suffices to show that in a 
small neighborhood of any C, S S'^ , the subspaces e_^p^fc(C), k = u, {v, q, a) (resp. ker r*^^(C), 
k = u, [v, q, a)) are close when — i^l + |^ — u| + |a| is small enough. Recall from (I3.17p (b) and 
(jl.47p that those spaces depend on the profile only through w{d). When — + — u| + |a| 
is small, we have w{d) ^ u, so 

(3.83) 4,^,(0 « 4,(0- 

3. Bounded frequencies. By compactness it suffices to show that for Q near some 
fixed Q, the corresponding spaces appearing in (13.8ip (b) are close for — i^l + |^ — u| + |a| 
sufficiently small. This is true for the spaces kerrfc((^), k = {u,ii), {v,q,a), since they 
depend on the profile only through w{d). The treatment of lE^(C) requires the conjugator 
P{z, C,, k) of Lemma l3.H where now we write k instead of p for parameters.. 

For / let IE^(C, k) be the stable subspace of ^(cxd, Q, k), which depends on the profile 
only through w{oo), and recall from the discussion below (|3.6p that 

(3.84) E^(C,A:) =P(0,C,fc)E-(C,A:), k = {v,q,a), {u,v). 
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The dependence of P{0,C,k) on the profile is not through w{0) alone but on the entire 
profile. However, recalling from Lemma |3. II and (the proof of) Corollary 13.31 that P{0,(,k) 
and E^(^,A;) depend continuously on the parameter k, we conclude from (j3.84p and the 
definition of the Evans function that 

for — idl + |(7 — ii| + \a\ sufficiently small. For C near 0, we replace lEi^(C) ^) by ]E(^(C) Pi k) 
(recall (jl.Sip ) in this argument. (This bounded frequency argument is the same as the proof 
of Corollary 13.31 but with the relevant parameters now explicitly identified.) □ 

4 Uniform Evans stability of small-amplitude layers for symmetric— 
dissipative systems 

In this section we prove Corollary 11.291 which shows that the uniform Evans condition 
Definition 13.4( b) holds for small-amplitude layers for symmetric-dissipative systems under 
several types of boundary conditions. By Theorem 11.281 it suffices to show stability of 
constant layers for symmetric-dissipative systems. 

4.1 The strictly parabolic case 

For clarity, we first carry out the simpler strictly parabolic case. 

Proof of Corollary \1.29\ in the case N = N' . 

Instability for rankTs > follows again by Proposition 12.51 combined with Lemma 
[031. 

1. Dirichlet boundary conditions. By Theorem 11.281 it is sufficient to prove 
stability of the constant layer w{z) = p. The matrix Y2{p) is now an invertible N x N 
matrix, so the boundary condition Y2{p)u{0) = h can be written u{0) = Y2~^{p)h := g. 
With A = ir+'j we consider decaying solutions u(xrf. A, rj) of the Fourier-Laplace transformed 
problem with coefficients evaluated at p: 

(a) XAqu + Adu' + i ^ Akr]kU 

(^•1) - Bddu" - i '^{Bdk + Bkd)r]ku' + ^ r]j7]kBjkU = 0, 

k^d j,k^d 

(b) n(0) = g, 

where the Aj are symmetric, Aq is positive definite, and Bjk is dissipative: 

(4.2) ftY^^Jk^J^k > 0\^\\ e > 0, for all ^ G M'^'. 

jk 
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Note that by the coordinate change 

(4.3) u = {Ao)-^^^w 

we can take Ao = I without loss of generahty, another advantage of constant-coefficients. 
We do this in the remainder of section HI 

By Remark 13.51 we must estabhsh the trace estimate 

(4.4) |n'(0)| <CA|<7|, 

where A ~ |1, r, 7!^*^^ + |t7|. In this case it is easy to treat frequencies of all sizes, so here we 
do not make use of the reduction to bounded frequencies effected by Corollary 13. 8i 
Taking the real part of the inner product of u with (j4.ip , we obtain 

3f?A(n, u)-^g- A^g - ^fig ■ Bddu'{0) + {u', ^B^du') 
^ - '^{r]ku, ?R.i{Bdk + Bkd)u') + ^ VkVjiu, ^Bjku) = 0. 

k^l j,k=id 

By (j4.2p . we obtain, extending to the whole line byu = Oonx< — ^ and u = {x\ri\ + l)u(O) 
on — ^ < X < 0, taking the Fourier transform, and accounting for errors introduced by 
extension, the Carding inequality 

{u, ^Bddu) + ^{r]ku, ^iiBdk + Bkd)u') 

(4.6) 

+ ^ %r?,-(n,3f?i?,fcn) > e{\u'\l + \r]\Ml) " CH\um\ 
j,k^d 

where 6 > 0. Here C|r/||ti(0)p is the error due to extension in the Carding inequality. It 
is an upper bound for the left side of (j4.6p . computed using the explicit formula for the 
extension of u and inner products in x < 0. Combining (j4.6p with (j4.5p . we obtain 

(4.7) (5RA + M')\u\l + \u'\l < C{\g\\u'm + (1 + |r?|)|5p). 

Similarly, taking the real part of the inner product of —u" with ()4.ip . we obtain after 
integration by parts, 

(4.8) m + M')\u'\l + lu'll < C ((|A| + |r?| + \r^\') \g\\u'm + (1 + h|)|tx'(0)|2). 
Using the Sobolev bound 

(4.9) |u'(0)|2 < |n'|2|u"|2 < CsA\u'\l + 6\u"\l/A, 

and the estimates on |n'|| and |n"|| coming from (14. 7|) and (14. Sp respectively, we obtain 

(4.10) A|n'(0)p<C(A2|5||^'(0)|+A3|5|2). 
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Since 

(4.11) A^\g\\u'm<CsA'\9?+SA\u'm\ 
we find that (j4.10p implies 

(4.12) A|m'(0)P < CA^Igp. 

2. Neumann boundary conditions. In place of (|4.ip (b) we now have the boundary 
condition 

(4.13) u'(0) = h, 

and we will continue to write u{0) = g m what follows. By Remark l3.5l it suffices to establish 
the trace estimate 

(4.14) A|m(0)| < C\h\ 

for a constant C independent of C,. To shorten the argument we now make use of the 
reduction to bounded frequencies |C| < -R provided by Proposition 13.81 In this regime the 
estimate (j4.14p is equivalent to 

(4.15) \u{0)\ < CA'^\h\, 

so we proceed now to prove the latter estimate. 
Again we use the Sobolev inequality 

(4.16) |u(0)p < \u\2\u'\l < 6\u\l + Cs\u'\l. 
Letting A' := 7^ + |?7|, for A' > c > we have from (fiTTll and (fi^ 

, , W\l < -^\9\\h\ + ^\g\' 

(4.17) 

\u'\l<CA'\g\\h\ + ^\h\'. 

Substituting into (|4.16p and taking 6 small enough, we absorb the terms involving g to 
obtain 

(4.18) |n(0)|2 < CA>|2. 

By assumption N'^ = N' = N and the fact that A^-|- = Nh — N'^ = 0, the boundary term 
—A^g ■ g in (j4.5p has favorable sign. Instead of (j4.7p we find 

(4.19) (3f?A + H^)\u\l + \u'\l + \um^ < Ci\g\\h\ + \r^\\g\^). 

For A' < c with c > small enough, we can absorb the terms involving g in (j4.19p to deduce 
()4.18p . This completes the proof of (j4.15p for bounded frequencies. 

3. Instability. The fact that the Evans condition fails for rankTa > N'i follows by 
Corollary 12.51 combined with Lemma ll.23[ □ 
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Remark 4.1. The argument of \GG^ establishing stability of small- amplitude Dirichlet pro- 
files in the non-constant coefficient, Laplacian-viscosity case required the weighted Poincare 
estimate 

r+oo 

(4-20) / \w'{z)\\u{z)\^ dz <\\zw'\\Lm)^+^)\\u'\\L2(^o^+oo), 

Jo 

established there for n(0) = 0. This estimate was also used in the one-dimensional treat- 
ments of JGS[ \R3\l . In the multidimensional case of general Bjk, the approach of jGO^ 
also requires a careful estimate of the error due to extension in the application of Carding 's 
inequality. These technicalities disappear in the constant- coefficient limit, a major simplifi- 
cation. However, it may be possible to use the argument of IGG] in certain cases involving 
variable multiplicities not covered by our structural assumptions. 

4.2 The partially parabolic case 

We now treat the case that Bjk is only semidefinite. 

Proof of CorollarylTM (N > N'). Instabihty for rankTa > N'l follows again by Corollary 
12.51 combined with Lemma 11.231 

1. Dirichlet boundary conditions. It is sufficient to consider the constant- 
coefficient equation 

(4.21) Xu + Adu' + i ^ AkTikU - B^u" - i^{Bdk + Bkd)Vku' + ^ r]jr]kBjkU = 0, 

ky^d kj^d j,kytd, 

where Aj are symmetric, Aq positive definite, and Bjk is now block-diagonal and dissipative: 

(4.22) 3? J2 ^j^kBfk > for C G \ {0}. 

For later reference we record the first and second components of (j4.2ip : 

(a) (ir + 7)-^! + i A]}r]kUi + i ^ Afr]kU2 + ^"n'l + Afu'^ = 
ki^d ki=d 

(4.23) + ^)''2 + i Afr]kui + i ^ Afr]kU2 + Afu[ + Af4 

k^d kj^d 

- Bfdu'i - iYiBfk + Blj)vkn'2 + Y WBfkU2 = 0. 

k^d j,ky^d 

By Proposition 13.81 we need only consider bounded frequencies |C| < -R; this greatly sim- 
plifies the analysis. 

Case (i) Totally outgoing flow. We first consider the totally outgoing case A]^ < 0, for 
which the boundary conditions in (|4.21[) are 

(4.24) U2{0) = g. 
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By Remark 13.51 a trace estimate that is equivalent to the (bounded-frequency) Evans hy- 
pothesis is 

(4.25) |ni(0)| + |n^(0)| <CA|g|, 

where A ~ |l,r, 7]^/^ -|- |?7|. We proceed to prove (j4.25p assuming (j4.22p and A^^ < 0. 
Pairing (j4.2ip with u we obtain the usual Friedrichs estimate 



-f\u\^ + \vf\u2\^ + W'2\^ + \uim^ < 

c{\u'2m\Mo)\ + \u2m' + \v\\Mo)\'), 



where | • | denotes an L?'{x) norm for interior terms and a C" norm for boundary terms, 
and we have dropped the subscript "2" on interior norms. Here the last term on the right 
represents the Garding error (from extension) as well as a boundary term from integration 
by parts. The |tii(0)p on the left is there because of the favorable sign oiAf. From (061) 
we obtain 

(4.27) |ni(0)|2 <C7A|<7|' + 5|n'2(0)|2. 

Similarly, differentiating (I4.2ip and pairing with u' we obtain 

7k? + |r/lV2l' + k2? + l<(0)l'< 
C7(|n'2'(0)|K2(0)| + (l + |r?|)K2(0)n, 

where again the terms on the right represent either Garding error or boundary terms from 
integration by parts. 

We now examine |ti2(0)|. First we have 

(4.29) l'"2(0)|^ < \U2\\U2\ < CsA\u2\'^ + jW^l"^. 
^Prom (j4.26p we find easily 

(4.30) C^AIu'al^ < CA^\g\^ + <5|u2(0)|^. 
We claim that the last term in ()4.29p satisfies 

(4.31) < CA2|5|2 + 5|4(0)|2 + 6\u,m''- 

With (|i:27P this wih complete the proof of (g^S]). 

To analyze the last term in (I4.29P we first use ()4.23p (b) to estimate |tt2(0)|: 

(4.32) |<(0)| < C {\X\\u2m + \7j\\u{0)\ + \u[{0)\ + \u'^{0)\ + |ryp|u2(0)| + |?7||n^(0)|) , 
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and then substitute in (|4.28p to get 



(4.33) 



C((|A||n2(0)| + |r?||n(0)| + |<(0)| + |n'2( 

+ \v?\uM\ + kll4(0)|)|n'2(0)| + (1 + mu'M?) < 
C(|A||U2(0)| + |r?||n(0)| + |??|>2(0)| + (1 + |r?|)|n^(0)|) 14(0)1. 



Here we have used the |n'^(0)p on the left to absorb a term, and then enlarged C. 
^From (I4.33P we have 

(4.34) < I ((|A||n2(0)| + |r?||..(0)| + |r?p|n2(0)|)|4(0)|) + |(1 + mu'MW 
and the second term on the right can be absorbed in (j4.29p . We have 

(4.35) ||A||n2(0)||4(0)| < <5A|^Z2(0)||u^(0)| < M^lyp + 6\u'^m\ 

The other terms in the estimate of ;^|ti2p are similar or easier to handle, so this concludes 
the proof of g^S]). 

Case (ii) Totally incoming flow. It remains to treat the totally incoming case > 0, 
with full Dirichlet boundary conditions u{0) = (lii(O), U2(0)) = g. A trace estimate that is 
equivalent to the (bounded-frequency) Evans hypothesis is 

(4.36) 14(0)1 <C7Ai|g|, 

where A ~ |l,r, 7]-^/^ + |??|. 

Making the same energy estimates as in the totally outgoing case, we find that the only 
differences are that (i) there now appears a term C|ni(0)p in the righthand side of (|4.26p . 
and (ii) there now appears a term C|4(0)p in the righthand side of (|4.33p . Difference (i) is 
harmless, since C|tti(0)p < C\g\'^ is of the same order as C|n2(0)p terms already appearing 
on the righthand side of ()4.26p . 

Difference (ii) can be handled by estimating |tix(0)| using (j4.23p (a) as 

14(0)1 < \{A\Y'\{C\C\\um + C7|4(0)|) < C((|A| + M)\g\ + 14(0)1) 

to see that C|4(0)P contributes terms of order C((|A| + If/D^I^P + I4(0)P)- Following the 
previous argument and using (|A| + \r]\f < A^, we obtain (14361) as claimed. 

2. Neumann boundary conditions. We now assume A'^^ = N' = rankTs. Consider 
first the totally outgoing case, so < and A^-|_ = = 0. The boundary conditions for 
()4.2ip are now 

(4.37) 4(0) =/i 
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and for bounded frequencies it suffices to show 



(4.38) |ui(0)| + |n2(0)| < CA2|/i|. 

Writing ^2(0) = g and A' = 72 + |t/|, since = we liave now in place of (|4.26p 

(4.39) 7|np + \v\>2\' + + H0)\^ < CHu'.miu^m + Hlu^m^). 

As before for A' < c with c > small enough, we easily absorb the terms involving g in 
(09I1 to obtain (08|) . 

To estimate n2(0) we use 

(4.40) \u2m^ <d\u2\^ + Cs\u'^\^. 
For A' > c> we have from (143911 and (14331) 

\n2\' < -^\h\\g\ + ^\gf 
(4-41) < ^1/^1 (|A|bl + \V\H0)\ + M'\9\ + (1 + |r?|)|/i|) 

<CA2|5||/^| + ^K0)||/.|+CM|/i| + ^|/.|^ 
^Prom (j4.40p . (j4.4ip we obtain, after absorbing some terms from the right 

(4.42) |n2(0)|2 < CA>|2 + C7|ni(0)||/i|. 

Substituting the estimate on |ui(0)| from ()4.39p into (j4.42p . we easily obtain (j4.38p after 
adding the estimates for ui{0) and U2{0) and absorbing terms. 

Consider finally the totally incoming case where A^^ > 0. We have rankTa = A'^^ = N', 
so N-^- = N\_ = N — N' and the boundary conditions for (I4.2ip are 

(4.43) ni(0) =51,4(0) =/i- 

Evans stability can fail in this case, even for small amplitude profiles. See Example 14.31 

□ 



Remarks 4.2. 1.) We note that the argument in the partially parabolic case is even simpler 
than the treatment of the one- dimensional case in IR3^ . thanks mainly to the reduction to 
finite frequencies and to trace rather than interior estimates. In particular, we do not require 
the "Kawashima-type" estimate used in IR3] to obtain an interior estimate on \u'i\. Nor do 
we require a weighted Poincare estimate of the type (j4.20p used in JR^. 

2.) In the small- amplitude analysis we expect that one may drop the constant-sign as- 
sumption on , substituting as in \R3j the assumption that that Fi be maximally dissipa- 
tive with respect to , i.e., that A^^ be negative definite on kerFi. For bounded frequencies 
the above arguments go through in this case essentially unchanged. For high frequencies the 
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exponential weights of IGMWZ6\. \Z3\l are not necessary in the small- amplitude case. It 
may be possible to modify the high-frequency argument of JGMWZBf to work without the 
constant-sign assumption. 

3.) In the case of Dirichlet boundary conditions (rankTs = Oj, the arguments of this 
section may be simplified still further by making use of Proposition below. This Proposi- 
tion together with Lemma \1.23\ yields uniform Evans stability for small frequencies |C| < r, 
r > sufficiently small. Using also the reduction to bounded frequencies ( Cor. I3.<gj] . it fol- 
lows that in order to show uniform Evans stability we must only show nonvanishing of the 
Evans function -D(C); that is, nonexistence of nontrivial solutions of the eigenvalue equation 
(j4.2ip with homogeneous forcing and boundary data, f = g = 0. Under these conditions the 
required estimates become almost trivial. 

Example 4.3 (A counterexample). Finally, we show that stability may fail in general 
for Neumann boundary conditions with A^^ = N' = rankTs, in the totally incoming case 
AY > 0. Consider the linear constant-coefficient system 

on X2 > with boundary conditions 

^tl 1x2=0 = 0) dx2U2\x2=0 = 0) 

under the assumptions 

6>0, b-a^<0. 

Then, > 0, and A^ = 1 = rankTs. Seeking layer profiles for this linear system, we have 
immediately, since these satisfy a first-order ODE in U2 with initial value at X2 = by the 
Neumann condition an equilibrium value, that these are exactly the constant solutions, from 
which we deduce immediately that the residual hyperbolic boundary condition is exactly 
ui = 0. 

Applying now Lemma fl.231 we find that a necessary condition for uniform low-frequency 
stability is satisfaction of the uniform Lopatinski condition for 

ut + AiUx^ + A2UX2 = 
with boundary condition ui\x2=o = 0, or 

(4.44) n/O for r {A2y\j + iT + irjAi)^ 

for each r, 77 G R, 7 S M^, where E_ as usual denotes the limit of the stable subspace as 
7^0"*". But, (j4.44p is clearly violated for 

r = (0,l)^, r/ = 1,7 = 0, r = 6/a. 
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in which case — (A2)~^(7 + it + ir]Ai)r = fir for fi = i/a. Continuing eigenvector and 
eigenvalue r = r('j), /x = /u(7) while varying 7 in the positive direction, we obtain by 
standard matrix perturbation theory, noting that r^A2 by symmetry is an associated left 
eigenvector, that 



verifying that r G E_. 

Alternatively, we may recall from |MZ2j that 2x2 two-dimensional hyperbolic constant- 
coefficient systems with both outgoing and incoming characteristics and for which Ai and A2 
do not commute satisfy the uniform Lopatinski condition if and only if they are maximally 
dissipative, i.e., ^2 < on the kernel span{(0, 1)} of the residual boundary condition, or 
< 0, to make the same conclusion without calculation. 

Note that this example does not yield one-dimensional low-frequency instability, and 
in fact is one-dimensionally stable. For, taking the real part of the inner product of u 
against the associated eigenvalue equation Xu + A2U' = U2, we obtain SRAIujl + \u2\2 = 0, 
yielding U2 = 0. Substituting into the ui equation, we find by direct computation that 
ui = e^'*'^ni(0) = 0. This gives one-dimensional Evans stability for bounded frequencies; 
the result for high frequencies follows by Proposition 13.81 

4.3 Maximal dissipativity of residual hyperbolic boundary conditions 

Before presenting calculations for example systems, we digress slightly to complete the 
picture of qualitative behavior. Transversality and uniform Lopatinski condition follow 
for small-amplitude profiles of symmetric dissipative systems by Evans stability combined 
with Lemma ll.231 under mild structural conditions on multiplicity of characteristics. Here, 
we note that they may alternatively be verified directly, and without any assumptions 
on multiplicity of characteristics. Proposition 14.41 yields the additional information that 
residual boundary conditions for small-amplitude layers of symmetric dissipative systems 
satisfy not only the uniform Lopatinski condition, but also the stronger condition of maximal 
dissipativity: SA^ < on the kernel of the linearized hyperbolic boundary conditions T^es, 
where 5 > is the symmetric positive definite matrix symmetrizing Aj. 

This result was established first in |BRa] for the case of symmetric, strictly parabolic 
-Bjfc. It was established in Lemma 4.3.1, [Met 4] for the strictly parabolic (not necessarily 
symmetric) case; see also Lemma 7 [BSZ] . The argument is based on dissipative integral 
estimates of a similar flavor to those used to establish uniform Evans stability in Sectional 

Proposition 4.4. Consider the class of symmetrizable dissipative systems (Definition \1.3\) 
and the class of decoupled boundary conditions (|1.4[) that are full Dirichlet in the parabolic 
variable U2 (N" = and maximally dissipative in the hyperbolic variable ui: that is, 
{SAdY^ is negative definite on kerPi. Then sufficiently small- amplitude noncharacteristic 
boundary layers are transverse, and the associated residual boundary conditions are maxi- 
mally dissipative. 






r'^A2r 
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Proof. 1. By continuity both assertions reduce to the corresponding assertions for the 
hmiting constant layer, say wlz) = p = Transversahty was already proved in 

Corollary 12.51 which applies also for the case of variable multiplicities. As in section H] (see 
()4.3p ) we may without loss of generality take S = Aq = I and Ai, = Ad {or u = {0, . . . , 0, 1)). 
We must show that A^ < on the kernel of the residual boundary condition Tres for the 
linearized hyperbolic problem at p. 

Let C^^p denote the manifold of states q near p such that there exists a profile W{z, q) 
satisfying (|1.16p with 

(4.45) (51,52) :=(Ti(pi),T2(p2)) 

and W{z, g) ^ g as z — > oo. Let Cy^p be the space of q such that there is a solution w{z., q) 
of the linearized profile problem (jl.lQp . (|1.20p with w{z,q) — > g as z ^ oo. The entries of 
G+{i^) in ()1.2ip are now evaluated at w^z) = p. It is not hard to check (see |Met4] . Prop. 
5.5.5) that 

(4.46) TpCi/^p — ^i^,p ' — kerrj.gs. 
The linearized hyperbolic problem at p is 



(4.47) 



vt + ^Aj{p)djV = f 
^1x^=0 £ TpCi/p. 



Therefore, we must show 

(4.48) ^rf < on C^,p. 
2. Set Fi := T[{p^) and define 

(4.49) AA = {(n,0) GM^ inGkerTi}. 
In Lemma 14.61 below we show 



(4.50) 4,p = AAeE_(A^i5, 



where E_(M) denotes the stable subspace of M. Since M C keri?^;^^ and A^ < on M 
(recall A^ < on kerFi), the result thus follows by Lemma 14.51 below. □ 

Lemma 4.5 ([Zl]). Let A and B be N x N matrices, where A is symmetric and invertihle, 
and B is positive semidefinite, ^{B) > 0, satisfying in addition the block structure condition 
ker B = ker ?R.(B). Given any subspace Af on which A is negative definite, then A is negative 
definite also on the subspace K-{A~^ B) © (MnkeiB), where E_(M) denotes the stable 
subspace of M. 
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Proof. Suppose that xq / lies in the subspace E_(y4 ^B) © (AAn keri?), i.e., 

Xo = Xl + X2 

where xi G K^{A^^B), X2 G {M H kev B). Define x{t) by the ordinary differential equation 
x' = A~^Bx, x{0) = xq. Then x{t) — > X2 as t — > +00 and thus 

lim {x{t),Ax{t)) = {x2,Ax2) < 0, 

t— >+oo 

with equality only if X2 = 0. On the other hand, 

{x,Axy = 2^{A-'^Bx,Ax) = 2^{Bx,x) > 0, 

yielding {xq,Axq) < {x2,Ax2) < 0. Thus, {xo,Axq) < unless X2 = and {?R.{B)x,x) = 0, 
in particular, xi), so that xi G ker3f?(B) = ker(i?) C ker A~^B. But, this is 

impossible, since xi G E,_{A~^ B). □ 

It just remains to prove: 

Lemma 4.6. Let Cu,p be as in (|4.46p . the space of q such that there is a solution w{z,q) of 
the linearized profile problem (I1.19p . (jl.20p with w{z,q) q as z ^ 00. Then 



(4.51) C^^p=MeE^{A^'B, 



dd) 



for N as in (|4.49p . 

Proof. 1. With w{z) = p the constant layer, define as in (j2.5p the N' x N' matrix 

(4.52) G,{p) := {Bjjr' {Af - Af{AY)-'Af) (p). 
A short computation shows 

(4.53) E^{A-'BM) = {(^'^^''^2^'"') ■■ r' G E„(G,(p))}. 

2. Consider the linearized profile equation (jl.l9p at p with W = (101,11:2,1113). For any 
{q^,q'^) G M^, this equation is easily integrated to yield a solution with {wi{z),W2{z)) 
(q^ ,q^) as 2; ^ 00: 

Mz) = -{A^^)-^Afe^^'^^P\G,{p))-\^ + 

(4.54) W2{z) = e'^''^P'>{Gd{p))-^r^ + q^ 

W3{z) = e'^^^'P^r^, where g E_{Gd{p)). 

Setting T[{p^)wi{0) = and T^(p2)ii;2(0) = we find 

T;(pi)gi = T;(pi)(yiy)-i^i^(G,(j.))-V2 
^ = -(G.(^.))-V^. 
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This gives 



(4.56) = {{q\q^) : T[{p') {q' + {A'.'y'Afq^) =0,q'G E_(G,(p))}. 



Together with (I4.53p . this imphes the result. 



□ 



Remark 4.7. A theorem of Ranch 'Ra]! asserts the converse result that any maximally 
dissipative boundary condition may be realized as the residual boundary condition associated 
with some symmetric dissipative viscosity. See also the interesting recent investigations 
of Sueur [Su] in which he establishes that any (nonstrictly) dissipative boundary condition 
may be realized as the residual boundary condition associated with some (not necessarily 
symmetric) dissipative viscosity. 



5 The compressible Navier- Stokes and viscous MHD equa- 



In this section we present computations of C manifolds for some classical symmetric-dissipative 
systems with various boundary conditions, including standard inflow/outflow conditions for 
Navier-Stokes. For the small-amplitude case the results of sections 12.11 and H] apply and 
give much information about when profiles satisfy the uniform Evans stability condition or 
transversality, and when the reduced hyperbolic boundary conditions, which are expressed 
in terms of C manifolds, satisfy the maximal dissipativity or uniform Lopatinski conditions. 
In a few cases we will say something about such properties for large amplitude profiles. 

5.1 Isentropic Navier— Stokes equations 

We start with the simplest case of noncharacteristic boundary layers for the isentropic 
compressible Navier-Stokes equations. In this case we are able to give a detailed description 
of possible boundary-layer connections, including large amplitude layers, and the resulting 
residual boundary conditions. 

Computation of residual boundary conditions. Consider the isentropic Navier- 
Stokes equations 



on the half-space y > 0, where p is density, u and v are velocities in x and y directions, 
and P = P{p) is pressure, and p > \rj\ > are coefficients of first ("dynamic") and second 
viscosity. We assume a monotone pressure function 



tions 



(5.1) 





(5.3) 



(5.4) 



P'ip) > 0. 
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Seeking boundary-layer solutions {p, u, v){y), and setting z := y, we obtain profile equa- 
tions 

(pv)' = 

(5.5) {puv)' = fiu" 

{pvy + P{py = ev", 

where ""' denotes d/dz, and e := 2p + tj. Integrating f^^, we obtain 

pv{z) = Poo^oo := "T-oo 

(5.6) pu' = puv - Pootioo^^oo 

ev' = pv^ + Pip) - (poov'L + P{poo), 

where {Poo,Uqo,Voo) '■= ip,u,v){+oo) and m := pv denotes momentum in the normal (z) 
direction. Within the set of allowable boundary conditions at z = in our abstract frame- 
work, we consider a subset consisting of linear conditions 

(5.7) TU{0)=g = {g\g\0) 

that are pure Dirichlet (N" = 0) or else mixed Dirichlet-homogeneous Neumann boundary 
conditions (A^" = A^' = 2), and which includes some of the most commonly used conditions. 

With = {{poo,Uoo,Voo) ■■ Poo > 0}, U{z) = {p{z),u{z),v{z)), Uoo = {Poo,Uoo,Voo), 

and Uo = {po,uo,vo), let 

Cr,g = {Uoo G : there exists U{z) satisfying ()5.6p together with 
U{+oo) = Uoo and the boundary conditions ()5.7p at z = }. 

For arbitrary Uq with pQ > note that the constant profile 

U{z) = Uo 

determines an element Uoo = Uq £ Cr,g, where g = TUq. The goal here is to determine 
Cr,g in some (not necessarily small) neighborhood of Uoo = Uq, and to understand how Cr,g 
changes with g. We consider separately the outflow and inflow case. 

Outflow with Dirichlet conditions. We first consider the outflow case vq < 0, with 
Dirichlet boundary conditions. We have < now, so N\_ = 0, Ni, = 2 and thus we 
prescribe only the parabolic variables on the boundary: 

(5.9) TU{0)=g = iuo,vo). 

Noting that moo = "lo < since f o < and pQ > 0, and rewriting the second equation 
in (15.61) as 



(5.10) fiu' = niooiu - Uoo) , 
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we obtain solutions 

(5.11) u{z) = Uoo 1- (1- — )e ^ 

satisfying the conditions (on u{z)) in (j5.8p with no restrictions so far except po > 0, Uoo £ M. 
Rewrite the third equation in ()5.6p using the first equation to get 



= moot^ + - rriooVoo + -t 



(5.12) moo(v - foo) + C^(yO, Poo) 



^oo '"-oo 



V V, 



oo 



moo (f - Voo), 

Woo 



where 



(5.13) c(p,poo):= >0 

P- Poo 

by assumption (j5.4p . Note that if we set 



(5.14) Coo = ^/P'{Poo) 

we have c{poo,Poo) = Coo- 

With (no, vq) fixed, consider first a candidate state Uoo for membership in Cr,g satisfying 

(5.15) Voo < 0, Uoo + Coo > 0, so |Uoo| < Coo. 

For z large we have v ~ Voo, P ~ Poo and thus 

fKtR\ rnooC^{p,poo) 

(5.16) moo > 0- 

Woo 

This implies that (j5.12p has no nontrivial bounded solution with limit Voo as z ^ +oo. In 
this case, the manifold Cr,g is an — A^4. = 2 dimensional manifold defined by the single 
condition 

(5.17) Voo = Vq, 

which is the right number of conditions (A^+ = 1) for the Euler equations at a value Uoo 
satisfying (j5.15p . Indeed, this is the classical (specified-normal velocity) Euler boundary 
condition corresponding to the (specified-normal velocity, no-slip) Navier-Stokes boundary 
conditions {u,v){0) = {uo,vo). 
In the case that 

(5.18) foo + Coo < 0, so |foo| > Coo (A^+ = 0), 
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we obtain the opposite inequality 



(5.19) 



"loo C^(/0, Poo) 



< 0. 



m, 



■oo 



Woo 



Thus, in this case we have nontrivial solutions (j5.12p with v{oo) = Voq. Provided vq is close 
enough to Vcxi, Cr^g is a full neighborhood of Uoo, as expected from the general theory. That 
is, there are no incoming hyperbolic characteristics and no boundary conditions imposed 
on the Euler equations in this case (dimCr,^ = N — = 3). 

Uniform Evans stability of the above layers in the small amplitude case follows by 
Corollary [TJHb). 

Examining further, we see that connections to such states Uoo exist for boundary data 
Uq such that 

(5.20) vl <vo< Voo or v^o < vq < v^, 

where the vi are the nearest rest points of (j5.12p to Voo- This is because all states v in the 
closed interval between vq and v^o then satisfy 



In the typical (genuinely nonlinear) case that P"{p) > 0, there is a single such rest state 
v"^ > "Voo- The states v"^ and t>oo then correspond to endstates of stationary viscous shock 
waves on the whole line. Layers with v < vq are in this case "compressive", consisting 
of pieces of the viscous shock from v'^ to v^o, while layers with v > vq are "expansive" 
(decreasing p), analogous more to rarefactions; see \MN\ ICHN Z] for further discussion. 

Remarks 5.1. There are two transitions worth mentioning. One is when the number of 
hyperbolic characteristics changes; that is, at the boundary where inequality (I5.15P changes 
its sense. Note that this transition has to do with the outer, hyperbolic solution Uoo CL^d 
so we cannot deduce that such a transition occurs from knowledge of the Navier-Stokes 
boundary data Uq alone, but must know the solution of the hyperbolic equation. 

A second transition has to do with the "inner", boundary-layer structure, when Uq goes 
out of range of ()5.20p . with Uoo held fixed satisfying (I5.19p . Consider a solution of the 
Navier-Stokes system with boundary data that includes states Uq with some in and some 
out of range of (15.20p . In this case our description of the solution as boundary layer plus 
outer smooth solution fails, and it is a natural question to ask what happens instead. This 
question has been answered for the one- dimensional case in IMNf . The resolution is that in 
this case there is a more complicated structure at the boundary consisting of boundary layer 
plus shock or rarefaction waves incoming to the domain: that is, a (nonsmooth) boundary- 
Riemann solution. 

Outflow with Neumann conditions. We next consider the case of homogeneous 
Neumann conditions in the outflow case: 





Woo 



(5.22) 
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Again consider a candidate state Uoo satisfying Voo < for membership in Cr,g- The 
vanishing of u'{0) and v'{0) imphes that u{0) and v{0) are rest points of (jS.lOp and (j5.12p 
respectively. Thus u{z) = u^o, v{z) = v^o, so t/oo can only be the endstate of a constant 
layer, and 

(5.23) Cr,g = {Uoo eR^ : poo>0,voo < 0}. 

In the case ()5.15p we have N+ = = 1, so N" = 2 > Nl; thus, by Corollary [23] these 
layers are not transversal, and thus are not even low frequency Evans stable by Lemma 
11.231 Observe that the "correct" dimension for Cr,g in this case is — A'^+ = 2. 

In the case when Uoo satisfies (IS.lSp . we have N+ = 0, so N" = = 2 and Cr,g has 
the right dimension. As expected there are no Euler boundary conditions and by Corollary 
11.29( b). the constant layers are Evans stable. 

Inflow with Dirichlet conditions. In the case of inflow vq > 0, so = 1, A^-|_ = 3 

and we prescribe 

(5.24) TU{0) = {po,uo,vo)=g. 

Let Uoo be a candidate state for Cr,g with Voo > 0. The equation for the transverse velocity 
u decouples as (|5.10p ; however, due to the opposite sign moo > 0, this equation now has no 
nonconstant solutions converging to Uoo- Hence, the transverse velocity is specified as 

Continuing as before, we find two cases, according as foo — Cqo ^ 0, i.e., 

"foo ^ Coo. 

In the first case A+ = 3 and we find as before that there are no nontrivial solutions of 
()5.12p . whence Voo = vq, and poo = ruoo/voo = mo/vo = pQ. Thus, only the constant layer 
is possible, and the induced hyperbolic boundary conditions are full Dirichlet, 

Uoo = Uo, Cr,g = {Uo} 

in agreement with the classical Euler conditions. 

In the second case A_|_ = 2 and there exist nontrivial connections on the range specified 
by (j5.2ip . The induced hyperbolic boundary conditions defining Cr,g are 

""oo = uo, ruoo = mo, 

specifying transverse velocity and momentum, an interesting variation on ()5.17p in the 
outflow case. 

Uniform Evans stability of the above layers in the small amplitude case follows from 
Corollary [LlUb). 
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Inflow with mixed boundary conditions. Finally, we consider the case of mixed 
Dirichlet-homogeneous Neumann conditions 



(5.25) rC/(0) = (p(0), n'(0), ^;'(0)) = (po, 0,0) = g 

for the inflow case. Again, u{z) = Uoo and v{z) = Uoo, since u{0) and v{0) are rest points of 
(jS.lUp and (j5.12p respectively. The mass equation then implies p{z) = poo, so in particular 
Poo = Po- Therefore, a candidate state C/qo can belong to Cr,g only if it is the endpoint of a 
constant layer, and we have 

(5.26) Cr,g = {UooeM:^:Voc>0,Poc = Po}, dimCr.g = 2. 
Consider the cases 

(5.27) i'oo - Coo ^ 0, 

and recah that Nb = 3 = N+ + Nl. In these cases N+ = 3, = and = 2, = 1 
respectively, so in both cases N" = 2 > N'^. Corollary 12 . 51 implies that these constant layers 
are never transversal; hence they fail to satisfy even low frequency Evans stability. Observe 
that in both cases Cr,g fails to have the "correct" dimension N — N^. 

5.1.1 Maximal dissipativity /uniform Lopatinski condition. 

The question of uniform Evans stability for large amplitude layers remains a mostly open 
question. For the one-dimensional isentropic case with 7-law equation of state, it has been 
shown numerically for Dirichlet boundary conditions that noncharacteristic boundary layers 
are stable, independent of amplitude [CHNZj. For the full, nonisentropic case on the other 
hand, it has been shown for Dirichlet boundary conditions that, even in one dimension and 
for 7-law equation of state, instabilities may occur [SZj . Here we show that the residual 
boundary conditions determined in the previous subsection are maximally dissipative (and 
thus satisfy the uniform Lopatinski condition) for all amplitudes. When F is a full Dirichlet 
condition (A^" = 0), this conclusion follows without any computation, since the residual 
boundary conditions by the analysis of the previous subsection have form independent of 
the amplitude of the boundary layer, and for small amplitudes are known (Theorem 14. 4p to 
be maximally dissipative. We carry out the computations nonetheless, to show how they 
work out in this simplest case. See |SZj for one-dimensional calculations involving more 
general, not necessarily decoupled, boundary conditions. 

1. Dirichlet conditions, outflow, Voo + Cqo > > Voq. The hyperbolic problem is 
symmetric, with one incoming, and two outgoing characteristics. Rewriting the isentropic 
Euler equations in U := {p,u,v) coordinates and linearizing about {poo,Uoo,Voo), we obtain 
Ut + A^Ux + A^Uy, where 





( Voo 





Poo\ 







^^oo 


° 




{clo/poo 





Voo/ 
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while the residual boundary condition (linear, in this case) is 



TresU = g 



r = (0 1) 



hence the kernel of the linearized boundary condition is the set of vectors (a, b, 0). Applying 
the positive definite symmetrizer S = diag{c^//9^, 1, 1}, we obtain 



which is evidently negative definite on the kernel of r^es; hence T^es is maximally dissipative. 

2. Dirichlet conditions, outflow, Voo + Coo < 0. In this case is negative definite, 
so the residual boundary conditions are automatically maximally dissipative. 

3. Dirichlet conditions, inflow, Voc, — Cqo > 0. Again, this is a trivial case, in which 
all hyperbolic modes are incoming, and the residual boundary conditions are TresU = U. 
So kerFj-es = {0}, and thus Tres is maximally dissipative. 

4. Dirichlet conditions, inflow, Voo — Cqo < < Voo- In this final case, we have 
prescription of momentum moo and transverse velocity Uoo, and the kernel of the linearized 
boundary condition is spanned by w := {poo,0, —Voo)- Computing 



we find that the restriction of 5*^2 to the kernel is again negative definite, so that the 
residual boundary condition is indeed maximally dissipative in agreement with the abstract 
theory. 

Remark 5.2. 1. We remark on the contrast between specification of velocity vs. momentum 
in outgoing vs. incoming case. (Specifying p or v along with u in the incoming case wouldn't 
be dissipative in general.) 

2. Neumann or mixed conditions. The only case in section [5n\ involving Neumann 
conditions that was Evans stable was the outflow case when Voo + Cqo < 0. Again, SA2 is 
negative definite so the residual conditions are maximally dissipative. In the remaining three 
cases involving Neumann or mixed conditions, maximal dissipativity fails. For example, in 
the inflow case (voo > 0) Tres = (1,0,0), and so SA2 is positive definite on kerr^es. 

5.1.2 Transversality of large amplitude layers 

We now show how to verify that transversality holds in the three cases where we constructed 
large amplitude layers: 

(a) outfiow/Dirichlet/voo + Cqo > 0; 

(b) outflow/Dirichlet/voo + Coo < 0; 

(c) infiow/Dirichlet/uoo — Coo < 0- 



(5.28) 
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It is helpful first to reformulate the definition of transversality geometrically. Let w{z) 
be a possibly large amplitude layer profile converging to Qoo '■= {Poo,Uoo,Vqo)- Set 

(5.29) Wo := (t/;(0),u;2(0)), Woo = iqoo,0) 

and observe that Wq G lies in both the stable manifold of Woo, denoted \V^{Woo), and the 
center-stable manifold of Woo, denoted W^^{Woo)- The conditions defining transversality 
of w{z) in Definition 11.111 can be rephrased: 

(i) r : Two {W{Woo)) is injective, 

(ii) r : Two i'W'iWoo)) is surjective. 

These may be recognized as the conditions that T be full rank on the stable and center- 
stable manifolds, respectively, of Woo, which correspond by definition to the following geo- 
metric version of transversality. 

Lemma 5.3 (Geometric transversality conditions). The transversality conditions of Defi- 
nition [T. 1 1\ are equivalent to the conditions that 

(i') the level set {W : TiyV) = T(Wo)} rneets the stable manifold of Woo transversally 
in phase space W = {w,w'l). 

(ii') the level set {W : T(W) = T{Wo)} meets the center-stable manifold of Woo 
transversally in phase space W = {w,vu'^). 

Remark 5.4. The use of the word "transversal" in (i') is perhaps nonstandard; it is used 
to suggest a minimal intersection. We mean {W : T{W) = T(l^o)} H W^{Woo) = {Wq}; 
that is, the intersection is a single point. 

The stable and center-stable manifolds are easily parametrized in each of cases (a) , (b) , 
and (c) above, so conditions (i) and (ii) can be checked explicitly. For example, consider 
case (c), where Nb = + = 2 + 1 = 3, the viscous boundary condition is 

(5.30) TU{0) = {po,uo,vo) ■.= Uo 

with Uq fixed, w{0) = Uq, uj{oo) = Qoo is fixed, and the dimensions of the stable and center- 
stable manifolds above are 1 and 4 respectively. Recalling the induced hyperbolic boundary 
conditions 

(5.31) Uoo = Uo, Poo'^oo = Povo, 
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we see that those manifolds can be parametrized as follows: 



Poo ^oo 

a 

a 




: a£ A] 



(5.32) 



\f{a,Poo,Voo)/ 

/ - \ 

a 



c 
a 




■.aeAcR,{b,c,d)eBc M^} 



\f{a,b,d)J 



where i? is a neighborhood of goo ^ is a neighborhood of vq determined by (j5.20p . and the 
function / (whose form turns out not to be important for verifying transversality) can be 
read off from (j5.12p . Independent vectors spanning the tangent spaces to WgiWoo) (resp. 
y^^'^ (Woo)) at Wq may be computed by differentiating the formulas (|5.32p with respect to 
a (resp. a, b, c, d). With these explicit formulas the injectivity and surjectivity conditions 
above are immediately obvious. 

We summarize this discussion in the following 

Proposition 5.5. The isentropic Navier-Stokes large- amplitude layer profiles described in 
cases (a), (b), (c) at the beginning of this section are transversal. 

5.2 Full Navier-Stokes equations 

Next we consider the full (nonisentropic) Navier-Stokes equations 

(a) pt + {pu)a, + {pv)y = 0, 

(6) {pu)t + {pU^)x + {pUV)y +Px = (2/i + 7])UxX + PUyy + (// + r?)^^?/, 

(c) {pv)t + {puv)x + {pv'^)y+Py = /i^xx + (2/U + 7])Vyy + (p + r])Uyx, 

(5.33) (d) {pE)t + {upE)x + {vpE)y + {pu)^ + {pv)y = 

KT^-X + HTyy + (^{2p + r])uUx + PV{VX + Uy) + tJUVy^ + 

{2p + ri)vvy + pu{vx + Uy) + r]vux ) 

^ / y 

where p is density, u and v are velocities in the x and y directions, p is pressure, and e and 

2 2 

E = e+'^ + \ ^'^^ specific internal and total energy respectively. The constants p> \r]\'>Q 
and «; > are coefficients of first ("dynamic") and second viscosity and heat conductivity. 
Finally, T is the temperature and we assume that the internal energy e and the pressure p 
are known functions of density and temperature: 



P = PiP,T), 



e{p,T). 
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It is clear for full (nonisentropic) gas dynamics that the uniform Lopatinski condition 
does not hold for arbitrary amplitudes, even for the simplest case of a 7-law gas. The 
proof of one-dimensional {rj = 0) viscous instability in [SZj for compressive boundary-layers 
of such a gas H was based on showing that the Lopatinski determinant (a multiple of A) 
restricted to the positive real axis could change sign as parameters, including amplitude, 
were varied; in particular, the determinant could vanish, yielding hyperbolic instability. So 
we cannot hope to show that the residual boundary conditions independent of amplitude 
as before. A brief examination reveals, likewise, that the residual boundary conditions 
are rather complicated to describe: in short, there appears to be no hope of other than a 
local analysis near the limiting constant layer. We carry this out below, and compute the 
linearized boundary condition for the linearized hyperbolic problem at a constant layer. 

The profile equations are, setting v := l/j, + rj: 

(pv)' = 
(puv)' = pu" 

{pv'^y+p' = uv" 

{pv{e + ^^^)+pv)' = Kr' + ]^{pu^ + vvy 

Let us introduce the unknown m = pv. Writing the equations in terms of the unknowns 
(m, M, t>, T) and integrating once gives the following equations with p = P{v, T) and poo = 



(5.34) (a) m = m 



00 

"2-00 



(6) u' = -^{U - Uoo) 

I \ I n-oo / \ , P Poo 

(c) V = (f-foo)H 

V V 

[d) T = [e - 600) - -Z—{U - Uoo) 

K 2k 

^00/ N.2 I Poo / X 
^^(^^-''^oo) H [V-Voo) 

Ik k 

The equation for X = {u,v,T) reads X' = F{X;moo, Xoq) where X^q and moo are param- 
eters. 



Case of an outgoing flow. Let U = (m, u, v, T) = (m, X) and consider an outgoing 
flow, V < 0, so Nb = N' + Nl = 3 + = 3. We fix a state Uq G Ug for which vq < consider 
boundary conditions that are just Dirichlet conditions on X: 

(5.35) r;7(0) = X{0) = Xq. 

^ Specifically, for 7 > 2 and 2fi + r; > k, in the inflow case v^o > > Voc — Coo with full Dirichlet 
conditions imposed at the boundary, instability was shown for layers consisting of a sufficiently large portion 
of a sufficiently large-amplitude 1-shock. 
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As usual we define Cr,Xo as the set of states (rraoo,-^oo) such that there exists a solution 
{m,X) of the equations ()5.34p on [0,+oo[, such that X{0) = Xq and U{oo) = (ttIoc^oo)- 
Clearly, Uq S Cr,Xo5 ™d i^^ar Uq we can write Cr,Xo the union 

(5-36) Cr,Xo = UmooCmoo.Xo 

where Cm^,Xo is the set of X^o S such that there exists a solution to equations 15.341 (b-c) 
on [0, +oo[ satisfying X(0) = Xq and X{oo) = Xoo- Corollary 12.51 shows that the constant 
layer U{z) = Uq is transversal, and thus Cr,Xo is a manifold near Uq. By the argument 
of [Sej, Lemma 15.2.5 (or alternatively, by an argument similar to our proof of Lemma l4.6p . 
the tangent space Tu^Cr^Xo is given by 



(5.37) 



TuoCr,Xo =^x ^-{DxF{Xo;mo,Xo)). 



Hence we are lead to calculate the stable invariant subspace of DxF{Xq] mo, Xq). One 
finds 



/ m 



(5.38) 



DxF{X;m,X) 







\ 



V i(Poo + 



me' 



J 



Since mo < 0, the matrix DxF{XQ;mQ,XQ) has at least one negative eigenvalue which 
is mo/fj,. It has exactly two eigenvalues with negative real part (in fact real negative 
eigenvalues) if and only if 

i(m + PO IP^ 
(5.39) det < 0. 



i(Poo + me;) ^e' 



T 



In that case E,~ (^DxF{Xo;mo, Xq)^ has dimension 2 and TuoCr,Xo lias dimension 3. Let us 
call A_ the second negative eigenvalue of DxF{XQ;mo, Xq). Then, the tangent space to 
Cr,Xo at the point {mo,Xo) is defined by the equations: 



(5.40) 



{■m,u,v,f) £ T(„o,Xo)Cr,Xo 



(i(mo + P^)-A_))t- + iP^r = 0. 



We do not know whether these boundary conditions have already appeared in the theory of 
the Euler equations, or if they have a special physical meaning. As we already know from 
the general theory, they are maximally dissipative for the Euler equations. 



5.3 MHD equations 

The equations of isentropic magnetohydrodynamics (MHD) are 
dtp + div(pn) = 

(5.41) < dt{pu) + div{pu^u) + \7p + H x cuilH = epAu + e{p + ri)\7divu 

^ dtH + curl(i? X n) = e^Ai? 
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(5.42) divH = 0, 

where p G M represents density, u G fluid velocity, p = p{p) G M pressure, and G 
magnetic field, with viscosities /i > \r]\ > and magnetic resistivity o" > 0; for simplicity, 
we consider here the case a = 0. When H = (|5.4ip reduces to the equations of isentropic 
fluid dynamics. 

The equations (|5.4ip are not yet in a form that satisfies the assumptions of section 
II. 1[ For example, the noncharacteristic condition is violated for every state {p,u,H). The 
equations may be put in conservative form using identity 

(5.43) H X cuvlH = (l/2)div(|iJp/ - 2H^Hf'' + i^divi^ 
together with constraint (j5.42p to express the second equation as 

(5.44) dt{pu) + div(pu*'u) + Vp + (l/2)div(|/7p/ - 2H^Hf' = epAu + e{fi + r/)Vdivn. 

They may be put in symmetrizable (but no longer conservative) form by a further change, 
using identity 

(5.45) cml{H x u) = {dwu)H + (n • V)H - {dwH)u - {H ■ V)u 
together with constraint (j5.42p to express the third equation as 

(5.46) dtH + {diYu)H + {u ■ V)H - {H ■ V)n = ae^H. 

Forgetting the constraint equation, we get a 7 x 7 symmetric system that satisfies all our 
structural assumptions except for (H4). 

Remark 5.6. 1) Define 

(5.47) c2=p'(p)>o, v = H/^, b = \ixv\, e = e/iei, 

cj := i (^c^ + + v^(c2-|t>|2)2 + 462c2) , 

c2 := i (^c2 + \v\^ - V(c2-|7;|2)2 + 462c2) . 
The boundary = is noncharacteristic for the hyperbolic part when 

(5.48) U3 ^ {0, ±V3,±Cs{n),±Cf{n)}, 

where Cs{n) and Cf{n) are the slow and fast speeds computed in the normal direction n = 
(0,0,1). Lemma 8.2 of JGMWZ^ shows that if we assume in addition 

(5.49) < |t;| / c, liisl > l^s], 

then Hypothesis (H4' ) is satisfied. For the stability analysis, which requires the construction 
of K -families of symmetrizers ( JGMWZ^ . Definition 3.5), we must use (I5.44p . (j5.46p . For 
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the purposes of deriving profile equations, describing C manifolds, and computing linearized 
residual boundary conditions, the form (j5.4ip . (j5.42p is more convenient to use. 

2) As the divergence-free condition is preserved by the evolution of the equations, we are 
free to ignore it in establishing stability of solutions. Results on instability, however, must be 
examined to check whether associated unstable modes are true, divergence- free instabilities 
or only apparent, spurious instabilities. 

With 1/ := 2/x + r/ the profile equations are: 
(5.50) 





{puz)' 


= 


{puius - 




= Pu'i 


{pU2U3 - 


H2HS)' 


= PU2 


{pul+p+^{Hf + Hl 


-Hi))' 


= ^^'i 


{Hius - 


- uiHs)' 


= 




- U2Hs)' 


= 






= 



By choosing the set of new unknowns m := pu^, Oi := — U1H3), (3 := {H2U3 — U2H3) 

H3, Ul, U2, U3 and integrating once, the system reduces to a system on u = {ui,U2,U3) 
where M := (m, a, (3, i/3) G and Uqo = (wiooi 1^2005 ^^Soo) £ are parameters: 



(5.51) 



M = 


Moo 




1 


m 




N OLH3 . 1 


Ul = 


— 


[ui - 


nioo) ( 






p -"3 




m 




^ 


u'2 = 


— 


[U2 - 


U-loo) ( 


At 




p -"3 


u's = 


m 
— 
u 


[U3 - 


U3 


+ 


1 

2^ 


{Hi- 


-Hfoo) + -^{hI 



1 



U300 

m 



Pi 



U300 



H. 



P- 

-M, 

)) + 
1 , 



■Ul 
^^3 
U2 
U3 



Uic 



U300' 
^^2oo ^ 
U300' 



2oo 



2u 

"+"^^^ and H2 



Hi - HI^). 



In the last equation it is understood that Hi = ^±Hiiii jj ^ P+U2m _ ^ 
calculations of section [5^2} the system has the form u' = F(u; M, Uoo)- The Jacobian matrix 
is 

a e/p 
a d/p 
e/p d/u b 



DuF{u; M, u) 
where M := (m, a, /3, H3) and 



(5.52) 



flU3 



UU3 



U3 



d :-- 



H2H3 

M3 
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The eigenvalues of DuF{u; M, u) are 



/. .ON ^ ^ a + b I (a-by + dP 

5.53 a and A± := — - ± W . + , 

2 y 4 /iz/ 

and we have 

(5.54) A+ A_ = ab- (e^ + 

Case of an outgoing flow. Let U = {m,a, (3, H^ju) = {M,u) and fix a state Uq = 
(Mo,«o) G l^d which satisfies the conditions (|5.48p and (j5.49p and for which the flow is 
outgoing, Us < 0, so Ni, = N' + N\_ = 3 + = 3. We consider boundary conditions that are 
just Dirichlet conditions on u: 

(5.55) TU{0) = u{0) = Uq. 

As usual we define Cr,uo as the set of states (Moo,itoo) such that there exists a solution 
{M,u) of the equations (|5.5ip on [0, +oo[, such that n(0) = uq and C/(oo) = (Moo,iioo)- 
Clearly, Uq G Cr,uoi aiid near Uq we can write Cr,uo as the union 



(5-56) Cr,«o = UmooCmoo,mo' 



where Ca,/oo,«o °f ^ such that there exists a solution to the last three 

equations in (|5.5ip on [0,+oo[ satisfying m(0) = uq and u{oo) = u^o- Corollary 12.51 shows 
that the constant layer U{z) = Uq is transversal, and thus Cr,Xo is a manifold near Uq. 
By the argument of section 15.21 the tangent space Tu^Cy^uq is given by 

(5.57) Tu,Cv,uo = X E- {DuF{uq; Mq, Xq)) . 

We proceed to compute the stable subspace of J^q := DuF{uq; Mq, Xq). 

The condition ()5.49p implies that J^q has at least one negative eigenvalue, namely a, 



d e_ 

three negative eigenvalues depending on whether 



with corresponding eigenvector (— ^j^^O). The matrix Tq will have exactly one, two, or 



(5.58) A_ > 0, A_ < < A+, or A+ < 

respectively. From ()5.57p we see that the corresponding dimensions of Cr,Mo ^' 6' 
respectively. By Proposition 12.41 the dimension of Cr,uo is N — N^, where is the number 
of positive eigenvalues of ^3. The eigenvalues of A3 are 

(5.59) Ao = Us, X±s = U3± Cs{n), A±2 = ^^3 ± I's, A±/ = U3 ± c/(n). 



so we observe that the cases (I5.58P correspond to the cases 

(5.60) lusl < Cs{n), Cs{n) < \u3\ < Cf{n), cj{n) < jusl 
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respectively. For example, consider the case A_ < < A-|-, which by formula ()5.54p occm's 
if and only if 

+ cP 

(5.61) ah- ^ < 0. 

A short computation now shows that the stable subspace of J-q is spanned by 
{ — ,-,0 , -,-,A_-a }, 

^^■^^^ , , b-a I {a -by e2 + d2 

where A_ — a = \ / 1 

2 Y 4 /izy 

The linearized residual boundary condition is thus expressed by 

/e d g2\ 

(5.63) (M, u) G rc/,Cr,«o ^u - (-{a- A_), -(a - A_), — + — = 0, 

where a,b,e,d are evaluated at Uq- 

Remark 5.7. The hypotheses and conclusions of Theorem \1.3(A on Evans stability and 

the existence of small viscosity limits apply to small amplitude layers for all the physical 

examples considered in section 5 except 

a)isentropic NS/inflow (vq > 0)/mixed Dirichlet- Neumann boundary conditions (j5.25p 
bjisentropic NS/outflow (vq < 0)/voo + Cqo > 0/Neumann boundary conditions ()5.22p . 
Additional examples for the full Navier Stokes and viscous MHD equations where The- 

orem [TTsU applies can be deduced from Corollary 1 1 . 29\( b). 



A Construction of Approximate solutions 

In this appendix, we give the construction of approximate solutions, following an approach 
similar to that used in [GGl IGMWZ7j . A new feature here is that we obtain global ap- 
proximate solutions on a domain with compact closure and smooth boundary. The same 
construction works on unbounded domains whose boundary coincides with a half-space 
outside a compact set. 

We seek high-order approximate solutions to 

d d 

{a)Ce{u) := Ao{u)ut ^^j(n)9jM dj{Bjk{u)dku) = 0, 

(^•1) j=i j,k=i 

{b)T{u,dTU^,duu'^) = (91,52,0) on dO, 



71 



which converge to a given solution u^[t,x) of the inviscid hyperbohc problem: 

£o(n°) = on [-TcTo] X J7 
u^t, xo) G C(t, xo) for (t, xo) G [-Tq, Tq] x OQ, 

where C{t,XQ) is the endstate manifold defined in Assumption 11.12] (see also Prop. 
Using the cutoff xi^) ™d the normal coordinates {xq,z) in a collar neighborhood of dfl 
defined in section 11.41 we look for an approximate solution of the form 

Uait,x)= e^U^{t,x,^) + e^'+^u'''+\t,x), 

(A.3) o<i<^^ 

W{t, x, -) = x{x)V^ it, xo, -) + u^t, x). 
e e 

Here satisfies ()A.2p and is given by 

(A.4) xo, Z) = W{Z, t, xo, u^(t, xo)) - u\t, xq), 

for a profile W{Z,t,xo,u^{t,XQ)) as in Assumption 11.121 

The V^{Z, xo,t) are boundary layer profiles constructed to be exponentially decreasing 
to as Z — > iboo. For the moment we just assume enough regularity so that all the 
operations involved in the construction make sense. A precise statement is given in Prop. 

E21 

A.l Profile equations 

We substitute (|A.3P into (|A.ip and write the result as 

Ad 

(A.5) Ye^T^{t,x,Z)\z=^+e^'R^^^'it,x), 

-1 

where we separate J^^ into slow and fast parts 

(A.6) P{t,x,Z) = F^{t,x) + G^{t,xo,Z), 

and the decrease exponentially to as Z ^ itcxD. 

The interior profile equations are obtained by setting the F^,G^ equal to zero. In 
the following expressions for G^ (t, xq, Z), the functions u^t,x) and their derivatives are 
evaluated at (i,xo). With W = W{Z,t,xo,u^{t,xo)) set 

Ut, xo, Z, dz)v := A,(^,^)(W)vz + iduA,(W) ■ v)Wz- 
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the operator determined by the hnearizing the profile equations about W, and 

d 

(A.8) Cov := Ao{u^)vt + ^Aj{u^)djv. 

i=i 



We have 



F-^{t,x) = 



(A.9) , d 

G-\t,xo, Z) = A,{W)Wz - —{B,{W)Wz) 

aZ 



(A.IO) 



G\t, xo, Z) = L(t, XQ, Z, dz)U^ - Q\t, xo, Z), 



where decays exponentiahy as Z — > +00 and depends only on [u^, V°). For j > 1 we 
have 

F^(t,x) = Cou^ - P^-\t,x) 
(A.ll) . ' ^ -^1 

G^{t,xo,Z)=Ut,xo,Z,dz)W+^ -Q^{t,xo,Z), 

where decays exponentially as Z — > +00 and P\ depend only on {u^, V^) for k < j. 
Note that dependence on the cutoff x{x) occurs only in the . 

In writing out the boundary profile equations, we note first that the boundary condition 
(jA.ip fb) is equivalent for e > to 

(A.12) T{u,edTu\ed,u^) = (51,92,0). 

With W {t, x, Z) = {W'^,W''^) always evaluated at {t,XQ,0) and 

(A.13) r{U^){v,0,vl):= [V.iU^y 

\KNdzv\ 

the boundary profile equations at order take the form: 
(A.14) T(Z^o,0,az^/0'2) = (51,52,0) (order e°), 

(A.15) T'{U''){U\0,dzU^'^) = (0,0,ci,3(t,xo)) (order e^), 

(A. 16) r(Z^O)(Z^^ 0, dzW'^) = (c,- 1, c,, 2, Cj, 3) (order e^,j > 2), 

where the Cj_fc(t, xq) depend just on the and their first derivatives for p < j — 1. 
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A. 2 Solution of the profile equations 

The solution of the profile equations given below assumes transversality of W{Z,u^{t,xo)) 
and the uniform Lopatinski condition, as well as the existence of a ET-family of smooth 
inviscid symmetrizers. Recall from Lemma 11.231 that the first two conditions both follow 
from the low frequency uniform Evans condition. 

1. The interior equations = and = and the boundary equation (|A.14p are 
satisfied because of our assumptions about and W{Z,t, xo,u^{t, xq)). 

2. Construction of {L(^,u^). We construct the functions h(^{t, x, Z) and u^{t, x) from 
the equations = 0, = 0, and the boundary equation (jA.lSh . will be a sum of 
three parts 

U^(t,x,Z) = U^ +U} + Ul , where 

(A. 17) \ ' ' / a CI 

Ul{t, x, Z) = x) + Vf^{t, xq, Z), k = a,b, c. 
First use the exponential decay of to find an exponentially decaying solution V^^(t, xq, Z) 

to 

L(t, xo, Z, dz)V} = Q\t, xo, Z)on ±Z>0 
(A.18) ^ 

^Oas Z ^ +00, 

and define u^{t,x) = 0. This problem is easily solved after first conjugating to a constant 
coefficient ODE using the operators P defined in Lemma l3. II 

Next, for fixed as above, use part (ii) of the definition of transversality (Definition 
11.111) to see that we can solve for Ul{t, xq, Z) £ S satisfying 



L(t, Xo, Z, dz)Ul = on Z > 



Recalling the definition of 5 from Lemma ll.lOl we see that lA^ has limits as Z ^ oo. Define 

ul{t,XQ) := lim Ul{t,xo,Z), 

(A. 20) 

Vb{t,xo,Z) :=Ul{t,xo,Z) -ul{t,xo), 

and let ul(t,x) be any smooth extension of it^(i, xq) to [— To,To] x Q. 

Finally, for an appropriate choice of ul{t,xo) we need U^{t,xo, Z) to satisfy 

Ut,xo,Z,dz)U^ = 
(A.21) T'{U'>){UlO,dzUt^) = 

lim U^{t,XQ, Z) = ul{t,X()). 

2— > + 00 
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According to the characterization of TqC{t,XQ) given in Remark ll.151 this is possible if and 
only if ul{t,xo) G T^o^^ .j.(j)C(t, xq). Thus, we first solve for ul{t,x) satisfying the linearized 
inviscid problem 

(A.22) ' 

u^{t,xo) G T„o(i_^.g)C(t,xo). 

This problem requires an initial condition in order to be well-posed. The right side in 
the interior equation of ()A.22p is initially defined just for t £ [—Tq,Tq]. With a cutoff 
that is identically one in i > —To/2, we can modify the right side to be zero in t < —Tq + S, 
say. Requiring to be identically zero in t < — Tq + 6, we thereby obtain a problem for 
that is forward well-posed since satisfies the uniform Lopatinski condition. Thus, there 
exists a solution to ()A.22p on [— ^,To]. This allows us to obtain U^{t,XQ, Z)) satisfying 
(fX2T]) and to define 

(A.23) V}{t,xo,Z) :=U'^{t,xo,Z)-ul{t,xo). 

By construction the functions (l/(^,u^) satisfy the equations = 0, = 0, and the 
boundary conditions (jA.lSp . 

3. Contruction of {U^,u^), j > 2. In the same way, for j > 2 we use the equations 
G^^^ = 0, = 0, and the boundary conditions ()A.16P to determine the functions {U^ ,u^). 
The corrector g^+i^(^^+i ig chosen simply to solve away an 0{e^^^) error that remains in 
the boundary conditions after the construction oiU^^ . 

In the next Proposition we formulate a precise statement summarizing the construction 
of this section. The regularity assertions in the Proposition are justified as in [ GMWZ4] . 
Prop. 5.7. Regularity is expressed in terms of the following spaces: 

Definition A.l. 1. Let (resp. H^) be the standard Sobolev space on [—Tq,Tq] x Q 
(resp. [-To, To] x d^). 

2. Let H" he the set of functions V{t,XQ, Z) on [— To,ro] x dO, x M_|_ such that V E 
C°°(R+,i/^([-To,ro] X dn)) and satisfies 

(A.24) |a|F(t,xo,Z)|H| < Ck,se~^^^^ for all k 

for some 6 > 0. 

Proposition A. 2 (Approximate solutions). Assume (H1)-(H6) (with (H4') replacing (H4) 
in the symmetric- dissipative case). For given integers m>0 and M >1 let 

T d ~\~ 1 

(A.25) so> m + - + 2M + —-. 

Suppose that the inviscid solution as in (|A.2p satisfies the uniform Lopatinski condition 
and that the profiles W{Z,u^{t,xo)) are transversal. Assume vP £ H'^" and u'^Iqq G 
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Then one can construct Ua as in ()A.3P satisfying: 

T{ua,dTul,duul) = (51,52,0) on dQ. 

We have 

(A.27) u^{t,x) G H'°-^^, V^{t,xo,Z) e ^^""^J', 



and i? (t, x) satisfies 



(^28) |(5t,5,.„,65,)°i?^"|i2 <C„ /or |a| <m+ ^ 

(6) \{dud,„ed,rR^'\L^ < /or |a| < m. 

In a collar neighborhood of dO., dx^ denotes an arbitrary vector field tangent to dO, . Away 
from such a neighborhood, can be a completely arbitrary vector field. 

B The Tracking Lemma and construction of symmetrizers 

In this appendix, we discuss further the high-frequency analysis of Section [3.21 in particular 
completing the proof of Theorem 13.61 by a treatment of the remaining (much easier) elliptic 
case. Recall that Theorem 13. 6| together with its easy consequences Corollary 13.71 and 
Proposition 13.81 allows us to reduce to considering only a compact set of frequencies in the 
stability analysis of many kinds of layers. We applied this result in showing, for example, 
that Evans stability of small amplitude layers follows from Evans stability of the limiting 
constant layer. 

In the process we prove a useful and previously unremarked relation between the track- 
ing lemma of |ZH|, IMaZ31 IPZ| and the construction of high-frequency symmetrizers as in 
[MZTl 'GMWZ3i I(;MWZ41 [GMWZ6]. In particular, we state the result (used implicitly 
in [GMWZ6J) that existence of a fc- family of symmetrizers implies continuity of decaying 
subspaces in the high-frequency limit (tracking). A similar argument was used in jMZ3] to 
establish continuity of subspaces in the low- frequency limit. 

B.l Abstract setting 

Consider a generalized resolvent ODE 

(B.l) U' -g{z,p,e)U = F, 

where ' denotes dz and {p, e) comprise frequencies and model parameters, in the limit as 
e ^ 0. We are interested in the situation, as in the high-frequency regime, that C is "slowly 
varying" in the sense that G' is small in an appropriate sense compared to Q. 
The following proposition gives one version of this notion in a simple case. 
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Proposition B.l. (a) Assume the M x M matrix G{z,p,e) is a C°° function of 
(B.2) (z,p,e) G [0,oo) X P X (0,1] 

for some parameter set P, and that \Q\ > ci > for ci independent of {z,p,£). Suppose 

(B.3) ^ ^ as e — > uniformly with respect to {z,p), 

and suppose 

(B.4) ^ f-f^ > Co > for Co independent of {z,p,e). 

Here, gap(M) denotes the spectral gap of a matrix M , defined as the minimum absolute 
value of the real parts of the eigenvalues of M . There exists a conjugating transformation 
T{Q) taking Q to block- diagonal form 

(B.5) TOT-^ = , mi+ > r]{z,p,e) and 3f?M_ < -r?(z,p,e), 

where r] := gap(^) > cq\Q\. For fixed co > as in (IB.4P the matrix T{Q{z,p,e)) satisfies 

(B.6) (a) |r| < c, (6) |r-i| < C7, (c) |T,| < c\g-M\g\, 

uniformly with respect to {z,p,e). 

(b) Setting U = T^^V in (jB.ip . we obtain 

V'-g{z,p,e)V = F, 

7^ g ■=gp + Q{z,p,£), 

where 5{z,p,e) = 5 := C\dzg\/\g\, and so 

(B.8) 6/r]^0as£^0 uniformly with respect to {z,p). 

Proof, (a). Consider the matrix fC := g/\g\ which has a spectral gap uniformly bounded 
away from by assumption. It follows by standard matrix perturbation theory |Kat| that 
there exists a smooth transformation T(/C) taking /C to block-diagonal form 

(B.9) TICP-^ = ^ ^ , > Co and KiV_ < -cq 

and satisfying (|B.6p (a). (b). Clearly the same T satisfies (jB.Sp . Moreover, with obvious 
notation we have 

(B.io) |r,| < \dT/d)c\ \d)c/dg\ \g,\ < c\g,\/\g\. 
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(b). Defining V hy U = T"V, substituting in (IRB . and using (|R9]) yields (iRTl) 
with 

(B.ll) F = TF, Q = -T{T-^)^, 

so by ([M \Q\<C\gz\/\g\. Thus, ([EH]) follows from the assumption (IM . 

□ 

For example, in the application of Proposition IB . 1 1 to the "elliptic zone" in section [BATI 
we will have 

(B.12) e = \G\ ~ Id, \Qz\ ~ Id, and so J ~ 1, r? ~ KL 

Depending on the application, r/ may be bounded, go to zero, or go to infinity as e ^ 0, 
and likewise |^|. In some sense this is artificial, since the scale of Q (and thus of rj) may be 
changed arbitrarily by rescaling the independent variable z; however, it is convenient not 
to have to rescale. 

More generally, following [MaZ3j . we take the existence of a transformation to form 
()B.7p - (IB.8j) as defining the notion of a slowly- varying coefficient Q. The construction of 
such transformations must in general be done quite carefully by hand, and does not follow 
by a simple argument like that of Proposition lB.lt see for example the treatment in Section 
7, [GMWZH] ■ and here in the proof of Theorem 13.61 of the general case away from the 
elliptic zone. Proposition IB.ll suffices to treat strictly parabolic systems, as pointed out in 
PZ1IZH1 IZ3]. 

B.2 Tracking 

For systems (jB.ip as above, we have the following version of the "tracking lemma" of |MaZ3j . 
The lemma implies that in the modified coordinates of (IB.7p . under assumption ()B.8|) . the 
subspace of initial data at zq G M of decaying (resp. growing) solutions of the homogeneous 
equation approximately "tracks" the stable (resp. unstable) subspace of the principal part 
Qp = blockdiag{M_|_, M_} evaluated at zq in the sense that one subspace approaches the 
other uniformly as e — > 0. 

In view of ()B.8P we will often suppress the dependence of 5 and r/ on {z,p) in what 
follows and simply write (5(e), rj{e). 

Proposition B.2 ( |MaZ3] ). Consider an approximately diagonalized system (IB. 71) with 
F = satisfying bound ()B.8p . 

(i) For all < e < eo, there exist (unique) linear transformations ^l{z,p) and <I>|(z,f3), 
with C°° dependence on {z,p,e) for which the graphs 

{{Zi,<^'^{z,p)Zi)} and {{^l{z,p)Z2, Z2)} 

are invariant under the flow of (]B.7p . The graphs consist precisely of the initial data 
at z = zq of solutions of ()B.7p that are respectively exponentially growing and decaying. 
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Moreover, the functions satisfy 

|^>^| < C6{e)/r]{e) ioi all z. 

(ii)In particular, the subspace E^{p,e) of data at z = for which the solution of (]B.7p 
decays as z ^ +00 is given by the graph {(<I>f (0,p)w_, : € C'^*"^*^"} and converges 
ase^O to E- := {(0,t;_) : V- G C'^*'"^-}. 

Proof. This can be proved by a contraction mapping argument carried out on the 
"lifted" equations governing the flow of the conjugating matrices see Appendix C, 
|MaZ3j for details. Proposition lB.4l provides an alternative proof of part(ii), which is the only 
part we use in this paper, in the case that Q exponentially approaches a limit as z ^ +00 
for each fixed p, e (not necessarily uniformly), which holds always in our applications. A 
related proof based on energy estimates/invariant cones appears in [ZH l IZlj . 

B.3 Symmetrizers 

With the same initial preparations (i.e., reduction to form (IB.7P ). one may also obtain 
directly bounds on the inhomogeneous resolvent equation by the method of Kreiss sym- 
metrizers p< | IMZl) . Consider (jB.ip on [0,+oo), augmented with some specified boundary 
condition 

(B.13) T{p,e)U = G, where \r{p,e)\ < C, 

uniformly for p £ P, e £ (0, 1]. The corresponding boundary condition for the conjugated 
problem (|B.7p is then 

(B.14) t{p,e)V := T{p,£)T-\0,p,e)V = G. 

Defining E^{p, e) as the stable subspace of the principal part blockdiag{M_|_, M_} evaluated 
at 20 = (i.e., E^ := {{0,v) : v G £^dimM^^^ defining the "frozen-coefficients Evans 
function" 

(B.15) ^(p,e) := det(^_,kerf), 

we have the following result. 

Proposition B.3. Consider the problem (jB.ip . ()B.13P under the assumptions of Proposi- 
tion \B.^ Assume this problem is 'frozen- coefficients stable" in the sense that \D\ > cq > 
for e > sufficiently small, uniformly with respect to p £ P. Then for some G > that can 
be taken independent of p £ P and e sufficiently small, 

(B.16) V^||C/||l2(r+) + \Um < C(||F|U2(M+)/V^ + |G|). 
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Proof. 1. In view of the properties (|R5]) . (lR6ll of the conjugating transformation T, it 
suffices to prove the estimate ()B.16P for V satisfying the conjugated problem ()B.7p and the 
boundary condition (|B.14p . 

2. Defining Sk ■= blockdiag{A:Id, —Id}, A; > 1, we have for e > sufficiently small, 

{i)Sk = SI, \Sk\ < k] 

{iiWSkQp > vie); 

{iii)SkV ■ V > A;|n+t)p — iXI-vp, 

Here Il± are the projections onto the unstable and stable subspaces of Gp{0,p,£). Writing 
V = (f + , f _ ) we have 

(B.17) U+v = {v+,0), n_ = (0,^;_). 

Thus Sk is a A;-family of symmetrizers in the sense of [GMWZ6] . 

Taking the real part of the inner product of —SkV with (|B.7p and integrating by 
parts, we obtain 

{l/2)SkV{0) ■ V{0) + {V,^{Skg)V) = ^{-SkV,QV) + ^{-SkV,F) 

<\Smvf + \\V\\\\F\\). 

^Prom this we may deduce using (i)-(iii) and 6/r] ^ that, for e > sufficiently small, 

(B.18) m+Vm^ + ir]/2)\\Vf < k\\V\\\\F\\ + |n_y(0)|2. 

But I det(ii^_, ker r)| > cq implies that |ri;| > c\v\ for v G where c = c(co,|r|) > 0. 
Thus, 

|n_y(o)| < |fn_y(o)| < \tv{o)\ + \tu+v{o)\ < \g\ + |f||n+y(o)|. 

Combining with (IB.lSj) . taking k sufficiently large relative to |r|, and using Young's in- 
equality to bound < (r//4)||y||2 + (|A;| V?/)||^f , we obtain the result. □ 

The next Proposition (more precisely, its proof) formulates the new observation that 
the conclusion of part (ii) of Proposition IB.2I is a consequence of the estimate ()B.16P ; that 
is, information provided by the tracking lemma may be deduced as a consequence of the 
basic symmetrizer construction. 

Proposition B.4. Consider an M x M system (jB.ip satisfying the assumptions of Propo- 
sition IBT^ Let E_(j), e) denote the subspace of initial data at z = for which the solution 
of 

(B.19) U' -g{z,p,e)U = 

decays to as z ^ oo. Let E_(j), e) denote the stable subspace of Q{0,p,e). Assuming 
(B.20) dimE_(p, e) = dimE_(p, e) forp e P, e small, 

we have 

(B.21) E_ {p, e) E_ {p, e) as e ^ 0. 
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Remark B.5. Condition ()B.20p holds, for example, if G exponentially approaches a limit 
as z ^ +00 for each fixed p, e (not necessarily as z ^ +00, by the conjugation lemma. 
Lemma This is always the case for our applications here. 

Proof. 1. Let E-(p,e) and E-{p,£) be the spaces appearing in Proposition IB. 2[ They are 
the exact analogues of E_ and E_ for the conjugated system V' — QV = 0. From (|B.5P and 
the properties of the conjugator T we have 

E_(p,e) =r-i(0,p,e)i?_(p,e) 
[ti.zl) _ , 

W._{p,e)=T-\Q,p,e)E_{p,e), 

so it is equivalent to show 

(B.23) E_ {p, e) ^ E_ {p, e) = {{0,v.) : v. £ c"'- } as e ^ 0. 

2. The proof depends on the fact that the estimate ()B.16p holds for any boundary 
condition satisfying the hypotheses of Proposition IB.3[ Let W C C*^ be any subspace 
transverse to the fixed space E^{p,e) and such that dimM^"*- = dim E^{p,e). Defining Tw 
to be orthogonal projection onto W'^, we have 

(B.24) ker Tw = W, rank Tw = dim E^{p,e). 

Since | det(-E'_,kerrvy)| > cq > 0, the hypotheses of Proposition IB.3I are satisfied by the 
conjugated system (IB.7j) with the boundary condition Tw- Assuming that V is any decaying 
solution of V — QV = 0, we have V^(0) G E^{p,e). From the estimate (|B.16P for the 
conjugated problem with F = we deduce 

(B.25) |y(0)| < C\TwV{id)\ for e sufficiently small. 

By (IB.20p and Lemma 11.211 this implies 

(B.26) |det(£^_(p,e),kerrvy| > c> 0, 

where c depends on C and iFv^l = 1- Therefore, W is also transverse to E-{p,e) for e 
sufficiently small. Since we are free to make different choices of W tranverse to E{p, e), this 
can only be true if ()B.2ip holds. 

□ 

As an immediate consequence of Proposition IB. 31 and Proposition IB. 21 part (ii), we 
obtain 

Corollary B.6. Suppose (jB.ip satisfies the assumptions of Proposition IBT^ and let 
(B.27) B(p,e) := det(E_,kerr). 

Then (IB.ip . (1B.13P is "uniformly Evans stable" in the sense that |ID)| > cq > for e > 
sufficiently small if and only if ()B.ip . ()B.13P is "frozen coefficients" stable in the sense that 

(B.28) lD)(p,e) := det(E_,kerr) > ci > 

for e sufficiently small. In either case the estimate ()B.16P holds. 
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Remark B.7. In this presentation, we have subsumed large parts of the usual symmetrizer 
construction into the preparatory transformations to the approximately block- diagonal form 
()B.7p . One of the main points of this development is to demonstrate that tracking and 
symmetrizer estimates are essentially automatic once we can reduce a system to form (|B.7p - 
(jB.Sp . Moreover, all high-frequency estimates derived in JZH[ \Z3l GMWZ4 , 1 GMWZ6^ were 



either obtained originally in this way, or may be rephrased in this form. 

Remark B.8. Another consequence of Corollaru \B.6\ is that uniform high-frequency Evans 
stability implies the maximal estimate ()3.8p for Q ^ £. Indeed, the estimate ()B.16p is 
equivalent to (jS.Sp in this application. 



B.4 Application to High Frequencies: Proof of Theorem 13.61 

We conclude by illustrating through explicit computations the application of these methods 
to the high-frequency analysis of Section 13.21 

Recall from ()1.43p the linearized eigenvalue equation Lu = f, where 

(B.29) L = -B{z)dl + A{z, C)d, + Miz, () 

with coefficients given by 

I3{z) = Bm{w{z)) 

d-l 

A{z, C) = Ad{w{z)) - ^Vj {Bjd + Bdj) (wiz)) + Ed{z) 



(B.30) 



d-l 

M{zX) = {iT + -i)Ao{w{z)) + Y,iVj{Aj{w{z)) + Ej{z)) 

i=i 

d-l 

+ X] ^3'nkBj^k{w{z)) + Eq{z). 

j,k=l 

The Ek are functions independent of (" which involve derivatives of w and thus converge to 
at an exponential rate when z tends to infinity. Moreover, we note that 

(B.31) eI^ = 0, = for A; > 0. 

With (|1.2p . we also remark that A4^^ does not depend on r and 7. 

As in (|1.45p we rewrite the linearized problem as a first-order system 

(B.32) d,U-g{z,C)U = F, r(c)t/|,=o = G, 

(B.33) 







g'' 


g''\ 










Id 






g32 


g33) 
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where U = {u,d,u^) = {u\u'^,d,u'^) £ C^+^' and C = (7,^,??). 
Here, 

^ _(^ll)-lA^12^ g32 ^ (g22)-l(_421^12 ^ _^22)^ 
^ _(^11)-1^12^ ^33 ^ (^22)-l(^21^13 ^ ^22)^ 

Note that Q^'^, Q^^ and are first order (hnear or affine in C), that Q^^ is second 
order (at most quadratic in Q) and that Q^^ is of order zero (independent of Q). We denote 
by Q'^ their principal part (leading order part as polynomials). We note that 

(B.34) gf{z, C) = Gf{w{z), C) when (a, h) + (3, 1), 

with 

Gii(n,C) = -(Ay(u))-i((7 + .rX(n) + E-=^%-A}i(n)), 

GlHu) =-(A\\u))-'Af{u) 

Gf{u,C) = {B^Hu))-'EtkliWB]l{u)), 

GfiuX) = -iB'Hu))-^ E •=! ^VAB]U^) + B^iv)). 

The principal term of Q^'^ involves derivatives of the profile w. Denoting by p = limz^+oo w{z) 
w{oo) the end state of the profile w, we note that the end state of is 

d-l 

There are similar formulas using the matrices Aj and Bj^k of (|1.5p . 
B.4.1 The elliptic zone 

We now prove Theorem l3.6l in the easiest case, which is for frequencies C lying in the elliptic 
zone 

(B.35) £ := {(r, 7, r?) : 7 > 6\C\ and |r/| > 6\C\} with 6 >0, 

with I el sufficiently large. (Recall that this is the final remaining case not treated in Section 
[321) 

In this case, the factors (I + 7) and A((^) = (^t'^ + + \r]\^^^^^ appealing in scaling (j3.1ip 
are both of order |(^|, and so we may replace (13. lip by the simpler rescaling 

Jc(u^u^u3):=(|C|^nMC|5nMC|-^^^) :=U 
Jc(5^5^5'):=(lClV,|ClV,|C|-V):=G. 
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If we define 

(B.37) T'f{u\u\u^) := {Tiu\T2U^ , K^u^ + \C\-^KTivW 

then ()B.32p may be written equivalently as 

(B.38) d,U-g{z,C)U = F, rf(C){7 = G' 

where 

/ a" g'' m'\ ..11 ^12 

(B.39) ^= ICIId := ' 



with obvious definitions of V"'^. Note that Q is of order one, while "P^^ is of order zero. 
Thus 

(B.40) g(z,C)=ap(z, + 0(1), ap=(^p' pS)=0(lCl), 

where the principal part Op is in fact homogeneous degree one in |(^|, and depends on z only 
through the profile w{z) and not its derivatives. 

Consider now the Fourier-Laplace transform of the upper triangular system 



j j 



written as a first-order system: 

(B.42) d.U-GutiOXW = 0. 

With U as above observe that we can write (jB.42p equivalently as 

(B.43) d,u-gp{o,ou = o- 

Proposition B.9. Suppose lies in the elliptic zone. LetK^(Q denote the space of initial 
data at z = of decaying solutions of dzU — g{zX)U = 0, and let E-(C) be the stable 
subspace of Qp{0,(). Then 

(B.44) ]E_(C) ^ E_(C) as |C| ^ 00. 

Proof. With p := \C\~^C ™cl e := |C|^^ we can write (with slight abuse) 

(B.45) giz,C)=giz,p,e) andgpiz,C) = gpiz,p,e). 
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The system (jR38]) thus has the form ([RT]) . By Lemma 7.3(i) of |GMWZ6| ) the matrix 
Q{z,p,e) has a spectral gap ??(e) ~ |C| = p Since 

(B.46) |g(z,C)|~|C|and|g,(z,C)|~|Cl, 

Proposition IB. II imphes there is a conjug ator T satisfying (|R5l) . (lR6l) . Hence 6{e) ~ 1 by 
([B:T0]) . (IB1T]1 . and thus 5/r] ^0 as e ^ 0. 

By Remark IB. 51 we can now apply Proposition IB. 41 to conclude that E_ ((^) approaches 
the stable subspace of Q{0,p,e) as e ^ 0. The 0(1) term in (jB.40p introduces an 0(e) 
difference between the stable subspace of Q{0,p,e) and that of Gp{0,p,£) (there is an 0(e) 
difference between the corresponding conjugators T), so we conclude ()B.44p . □ 

Proof of Theorem \3.6\ for the elliptic zone. Here we use the notation of the previous Propo- 
sition. 

1. For C, & £ and as in (1B.38|1 define a slightly modified rescaled Evans function 
(B.47) Df(C) = det(E_(C),ker Tf). 

Using Lemma 11.211 in the same argument that showed the equivalence of estimate (j3.10p 
and the uniform Evans stability condition Definition 13. 4|, we see that for C ^ ^) 

there exist cq, R such that |-D'*''(C)I > co for |C| > i? 
there exist ci, i?' such that \Df{C,)\ > ci for |C| > R. 

2. Defining a frozen-coefficient Evans function 
(B.49) D'^'iC) = det(E_(C),kerr 



e )i 



we conclude from Corollary IB. 61 and Proposition IB. 91 that the conditions (|B.48P are equiv- 
alent to 

(B.50) there exist C2,R" such that ]d'''{C)\ > C2 for |C| > R" ■ 

Recall that IE_((^) is the stable subspace of Gp{0,C), which defines the problem (IB.43P 
equivalent to the upper triangular system (|B.4ip . 
3. To complete the proof we write 

(B.51) Tf = (ri,rf) 

and observe that 

^-n-(Y')H.-V 

^Prom the form of Qp we see that lE-(C) consists of vectors of the form 



(B.52) ( , va € SiGl') and f * ] , v, e S{Vf), 
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where S{M) denotes the stable subspace of M. After performing obvious column operations 
on D^'^{Q to express the determinant in upper block triangular form, we thus obtain (up 
to a sign) 

(B.53) D'^iC) = det (^(^^^1), ker Ti) x det (^(pf), ker Tf ) . 

Although r^^(C) here differs slightly from r^'^(C) in (|3.2Up . the argument that established 
the equivalence ()B.48P allows us to conclude that D^(C) (|3.20p is bounded away from for 
\C\ large if and only if the second factor on the right in ()B.53p is. The first factor equals 
D^iC) (|3-20p . so this completes the proof. 

□ 

Remark B.IO. For coupled boundary conditions, the above analysis shows that the high- 
frequency Evans condition is more complicated, involving at least the upper triangular system 
()B.4ip rather than the decoupled (jl.Sp . 

B.4.2 The general case 

In the remaining frequency regimes 

Cs:={C:0<l<m}U{C:\r]\<5\C\}, 

S > sufficiently small, the reduction to form ()B.7p is considerably more complicated, 
in particular involving multiplication of hyperbolic modes by an exponential weighting 
function to obtain definiteness of M±. Moreover, the resulting terms M± depend in hyper- 
bolic modes both on derivatives of the profile and on the chosen exponential weight. Thus, 
it is no longer true that the effective principal part Q after rescaling/appropriate coordinate 
transformations involves only frozen coefficients of the triangular system (IB.4ip as in the 
previous case. However, it is still true that the associated stable (resp. unstable) subspaces 
of Q depend only on those coefficients: indeed, only on the coefficients of the fully decoupled 
system (j3.17p . Thus, the above conclusions about the rescaled Evans function and uniform 
high-frequency stability (based on the form of the stable subspaces and not on M±) hold 
true in this case as well. The maximal estimates follow likewise from Proposition IB. 31 as 
before, once form ()B.7p has been achieved. For detailed treatments, including in particular 
an implicit reduction to form (IB.Tp . see Section [32] of this paper and Sections 7.2-7.4 of 
[GMWZGj . 
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